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ABSTRACT 


The objective of this thesis was to contribute to the development of a second-order 
approximation to the steady and oscillatory lifting flow past low aspect ratio wings by establishing 
the validity of a NASA-Ames developed incompressible three-dimensional flow panel code named 
PMARC, which would subsequently be used to test the range of validity of this second-order theory. 
The steady state configurations modelled in the validation process include swept back, delta and 
rectangular wings, an F5 wing and three equivalent bodies of revolution. Oscillatory cases were also 
run with one delta and the F5 wings and their respective equivalent bodies of revolution, as well as 
with one spindle. In steady state, comparison with experimental and theoretical data proved 
PMARC to be very accurate for lift and pressure calculations, but revealed a discrepancy in the 
velocity distribution calculation around delta wings. This finding was corroborated by applying the 
slender body/slender wing theory (Oswatitsch-Keune theory) to the delta wing. The unsteady state 
results are presented, but their validation is left for future work. As part of its primary objective, this 
thesis also presents a computer code that generates the F5 wing equivalent body of revolution from 
its chordwise section definition. In order to be used as base software to the second-order theory, 
PMARC will have to be corrected, or a new software will have to be validated. 
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I. INTRODUCTION 


Oyer the past quarter century, great advances have been made in the computation of the 
aerodynamic forces on aircraft and missiles. These advances were made possible by the 
development of high-speed computers which, in turn, made possible the numerical integration of 
the compressible flow equations. However, the numerical computation of the flow field around a 
typical military aircraft or missile configuration is still a very time consuming and expensive task 
because it requires the development of a complex three-dimensional computational grid and the 
solution of a system of nonlinear partial differential equations by means of a suitable numerical 
integration scheme. Furthermore, the design of aircraft or missiles requires not only the 
determination of the aerodynamic forces in steady flight but also of the transient and oscillatory 
forces in order to assure the flight vehicle's dynamic and aeroelastic stability and to predict the 
dynamic loads due to wind gusts. For these reasons there exists a continuing need to develop 
approaches which are sufficiently accurate, yet are computationally more efficient than full-fledged 
three-dimensional solutions of the inviscid compressible flow equations (Euler equations) or of the 
viscous compressible flow equations (Navier-StokeS equations). 

It is the objective of this thesis to contribute to the development of a new approach to 
determine the steady and unsteady airloads on missiles and on low-aspect ratio aircraft 
configurations typically used for high-performance fighter aircraft. This approach is based on an 
extension of the well-known slender-body/slender-wing theory developed by Oswatitsch, Keune, 
R.T. Jones and others in the 1940's and 1950's (see e.g. the book by Ashley and Landahl, Ref. [8]). 
Oswatitsch and Keune showed that, in a first approximation, the non-lifting flow over bodies of 
revolution and low aspect ratio configurations can be obtained by the superposition of two flows, 
namely the flow in the planes normal to the ffeestream direction (the cross-flow) and the "spatial 
influence" which accounts for the contribution upstream and, in subsonic flow, also downstream of 
the cross-flow plane under consideration. The cross-flow can be treated as a purely incompressible 
flow (even though the missile or aircraft flies at transonic or supersonic speeds), whereas the spatial 
influence is Mach number dependent but has the remarkable feature that the spatial influence is the 
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same for all configurations which have the same cross-sectional area distribution in the streamwise 
direction. The lifting flow, on the other hand, can be obtained by merely considering the cross-flow, 
hence reducing the analysis to a purely incompressible flow analysis for high-speed flight vehicles 
of low aspect ratio. 

M.F. Platzer has published a NASA technical memorandum (NASA TM X-53295, April, 
1965) on his development of a second-order approximation to the steady and oscillatory lifting flow 
past low aspect ratio wings which exhibits the same fundamental division into an incompressible 
cross-flow and a Mach-number dependent spatial Influence. To test the range of validity of this 
second-order theory it is necessary to compare it with an exact three-dimensional flow solution. 
Fortunately, such a solution has become available for low speed (incompressible) flow in recent 
years with the NASA-Ames developed PMARC code. Therefore, in this thesis, the following 
specific objectives were attempted, namely: 

1. Demonstrate and document the PMARC code, 

2. Compare the PMARC solutions with other available solutions in order to establish 
its validity, 

3. Compute the incompressible flow over a delta wing at zero angle of attack using the 
Oswatitsch-Keune theory, 

4. Obtain the equivalent area body of revolution of an F-5 wing, 

5. Compare the Oswatitsch-Keune theory with the PMARC predictions, and 

6. Generalize the non-lifting theory to incompressible flow past oscillating bodies of 
revolution and low-aspect ratio wings. 

These tasks will be described in the subsequent chapters. 
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II. SOFTWARE OVERVIEW 


A. THEORY BEHIND PMARC 

PMARC stands for Panel Method Ames Research Center. It is a potential flow panel code 
developed at NASA Ames Research Center to numerically predict flow fields around complex 
three-dimensional geometries. The general idea behind this panel code is to cover the body under 
study with panels over which there is an unknown distribution of doublets and sources. (Recall that 
the doublet and source flows are two of several elementary flows used to model the flow field 
around arbitrary shapes.) These unknovras are solved through a set of simultaneous linear algebraic 
equations generated by calculating the induced velocity at each panel's control point and applying 
the flow tangency condition. 

Two kinds of simplifications were made to allow this software to remain computationally 
cost effective. The first one has to do with the strength of the elementary flows (also called 
singularities). In PMARC, the strength of the singularities is constant over each panel allowing 
factorization of the doublet and source strengths in the thus greatly simplified analytical solution of 
the flow. The second kind of simplification is in the assumption of an inviscid, irrotational, and 
incompressible flow field, thus satisfying Laplace's equation: 

V 2$=0 ( 2 . 1 ) 

With these simplifications, PMARC proceeds to solve for the potential at any point P by 
applying Green's Theorem, which results in the following integral equation: 

dS - [(4)«(V®-V®^] dS (2.2) 

4nJJs^w*s„ f A-kj Js*w*s„ f 
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where r is the distance from the point P to the element dS on the surface and n is the unit vector 
normal to the surface pointing into the flow field of interest. In this equation the first integral 
represents the disturbance potential from a surface distribution of doublets with strength 'phi-phi,' 
per unit area and the second integral represents the contribution from a surface distribution of 
sources. As shown on Figure 2.1, phi corresponds to the external region (or flow field of interest), 
and phi, is the internal region (or fictitious flow). This arrangement is required to model external 
flow problems such as a wing in a uniform stream. 



(From Ref.[l]) 

PMARC actually solves a simpler form of Green's Theorem than the one presented above. 
The simplifications are made possible by noting that at the surface at infinity, the perturbation 
potential due to the configuration is essentially zero, leaving only the potential due to the uniform 
onset flow. Also, it is assumed that the wake is thin and there is no entrainment so the source term 
for the wake disappears and the jump in normal velocity across the wake is zero. 
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The theory described above makes for three major steps in the PMARC solution process. 
First, it solves for the doublet strengths leading to known values of the singularities on each of the 
panels. Second, it computes the velocity by differentiating the velocity potential expression in the 
appropriate direction for each tangential component of velocity. Using this information it calculates 
the resultant velocity at each panel control point. With this result it finally derives the pressure 
distribution using Equation (2.3). The forces and moments on the body are then obtained by 
performing a surface pressure integration. 



A? 


(2.3) 


B. PMARC OPERATING MODE 

Instead of giving all the details pertaining to the operation of the PMARC software within 
this section, only the generalities will be presented. For practicality and ease of understanding, the 
information and instructions for construction of a specific wing configuration will be presented in 
the software validation sections, each time a new configuration is introduced. For additional details, 
the reader should refer to Refs. [1], [2], and [3]. Installation instructions are also found in Ref. [1]. 

1. Generalities 

There exists more than one way to construct a structure in PMARC. One is to make use of 
three levels of coordinates, as shown in Figure 2.2, in order to join components built separately. 
Another way would be to build the whole structure at once referring all points to one coordinate 
system. 
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Figure 2.2. PMARC Hierarchy of Coordinate Systems. 
(From Ref. [2]) 


But no matter which method is preferred, the physical representation of each component 
is made with sheets of panels called patches used to solve the required simultaneous flow equations 
for that particular group of panels. The panel density of these patches and the manner in which all 
the patches are joined together are defined within the input file as will be shown in later sections. 
The shape of theses patches is given by providing a point definition, or path, around or along which 
the patches will be laid down or wrapped. Figure 2.3 shows how a patch is rolled up on top of itself 
to make a surface in which two of the four sides are joined in a seam. There are a few rules about 
the numbering of each panel, and although most of them are transparent to the users upon 
construction of a model, others shall be known for output interpretation purposes. These rules will 
be explained later. 
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Figure 2.3. PMARC Patch of Panels. 
(From Ref. [2]) 


The previous figure is a representation of hoAv wings are modelled. Again, PMARC offers 
some flexibility as there are two ways to define a wing section. The first one is done by entering a 
series of points that define the wing cross-sectional shape. PMARC will then wrap the patch on top 
of that definition. The second method will automatically generate the point coordinates of a NACA 
Four Digit Series airfoil. All that is required are the airfoil's maximum thickness ratio, and the 
maximum camber and its location. A later section will show how to use these two options. 

Once the shape of a section is defined, the third dimension of the component has to be 
modelled. Again, several options are available to this end. One is to position the leading edge origin 
of each end of the component, for instance the root and the tip cross-section of a wing, in some 
chosen coordinate system, and then to scale each of these ends to the desired size. Another method 
is to define the cross-section of each end of the component in relation to a coordinate system. This 
way, each cross-section is built originally to its desired size, instead of being scaled from the non- 
dimensionalized definition as with the previous method. Depending on the application, each of these 
methods, or a combination thereof, show significant advantages. 







PMARC incorporates routines to automatically generate wing tips, whether they need to be 
flat or rounded. Also, by toggling the value of certain variables, complete components can be copied. 
A multitude of other features are available among which is found the wake generation and the 
streamline generation routine. The interested reader should refer to Refs. [1], [2], and [3]. 

2. PMARC Input File 

Several decisions have to be made when building a geometry with PMARC. The number 
of panels spanwise and chordwise and their distribution are two important ones. More panels means 
better accuracy, but also means proportionately longer computational time. The optimal number of 
panels should be sought to make the results sufficiently accurate. The panels can be spaced 
following four types of distribution: full cosine, two types of half cosine, and equal spacing. These 
allow the user to concentrate more panels in areas of expected greater gradients. These decisions will 
be justified in the validation section each time a new configuration is tested. 

Other information, geometry specifics, and test conditions are entered into PMARC via its 
input file. It is called DATA5 and is FORTRAN based. Over 200 fields or variables are required 
to be filled for PMARC to successfully read the file. Although not all the fields are pertinent to all 
cases, an uncompleted field will cause a software crash. The following paragraphs will laconically 
describe the DATA5 input file, identifying which fields should be filled with correct values to obtain 
correct output in the case where we are interested in the Cp, Vx, Vy, Vz, V, CL, and Cl but without 
wake or streamline generation. The reader should refer to Figure 2.4 before continuing. Also, each 
field of the DATA5 input file is represented by an acronym. The definition of these acronyms along 
with the valid entries are listed in Ref. [1] starting on page 33. 
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Figure 2.4 PM ARC Input File Sample. 


&VS3 Xl= -0,1000, YU 1.5000, ZU 0.1000, MPTU 0, 

iiVSil X2= -0.1000, Y2= 1.5000, Z2» 0.1000, MPT2= 20, 













a. Basic Input Section 

This section allow the user to entitle the file, prescribe specific job requirements, 
control some solver and time-step parameters, define some of the testing conditions as far as the 
movement of the model and the ffeestream are concerned, and also to input reference dimensions. 
All the above make up for sixty fields. Not all of them require actual testing values in order to 
obtain valid output for the parameters listed in subsection 2 above. These fields can be filled with 
any number without altering the results. In the steady case, the critical fields are the type of run 
(LENRUN) that determines if a full flow calculation will be performed or just a geometry 
generation, the angle of attack (ALDEG), and the reference surface area (SREF), the latter being 
used in the local lift coefficient calculations. Other fields become critical when dealing with 
unsteady state cases. They will be discussed later in the oscillation chapter. Two of these fields 
(NTSTPS, DTSTEP) also require some attention in the steady state cases as will be seen. 

One should always study the output carefully as some of the results might be the 
consequence of fields that were filled with arbitrary values instead of values corresponding to the 
test case. A prudent habit for the analysis of any kind of wing or wing-body combination fixed in 
space is to fill the following list of fields with the information corresponding to the object under 
study: LENRUN, VINF, VSOUND, ALDEG, CBAR, SREF, and SSPAN. The other fields can 
remain the same as the values used in example files presented in the PMARC manual (Ref. [1]). 

b. Input Section for Surface Geometry 

This portion of the input file includes the headings &ASEM 1, &ASEM2, &COMP 1, 
&COMP2, &PATCH1, &SECT1, and &BPNODE. For the example analysis described, the first four 
sets of fields don't change from what is shown on Figure 2.4. They will change if two or more 
independent components are studied. The last three sets (&PATCH1, &SECT1 and &BPNODE) are 
critical in the construction of components and the definition of its panel layout. Each individual 
field within these sets has to be entered with values specifically corresponding to the model being 
built and the conditions of the tests. As will be seen later, the sets of fields listed above are 
duplicated as often as required to define more complex components or assemblies. Also, &PATCH2, 
&PATCH3, &SECT2, and &SECT3 are other sets of field required in certain cases. 
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It becomes a complex matter to identify which cases require which ones of the sets 
of fields discussed above. The uninitiated should study the various test cases in the validation 
sections to acquire a broader view of the many possible applications. 

c. Time-stepping Wakes and On-body Streamlines 

Wake analysis with PMARC is a challenging undertaking. If wakes are to be 
modelled, each of the fields in the wake portion of the input file is important. These fields lose all 
significance when the wake routine is turned off (IDWAK = 0). The same comments are true for the 
on-body streamlines. 

3. PMARC Output Files 

Once the PMARC has successfully read the input file DATA5, it performs the flow 
calculations and generates two output files named DATA6 and DATA22. The latter is an unreadable 
file that contains all the information required for the visualization of the various parameter gradients 
along the body under study. This visualization is made possible via a satellite software named GVS. 
General information on GVS can be found in Ref [3] and a detailed user guide is also available in 
Ref [18]. 

DATA6 is the output file formatted for interpretation by the user. Its length could range 
fi’om a few kilobytes for the simple geometry computation to several megabytes for full flow 
computation using several time-steps. DATA6 is a well formatted and structured output file. The 
following provides a succinct description of its structure. 

a. Input list for Basic Data 

The top of that file reiterate some of the important geometry and test conditions 
prescribed by the user within the input file, DATA5. It also lays out the position of the assembly and 
component coordinate systems. This portion of the file ends with the cross-sections definitions and 
a list of the basic patch data. 
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b. Time-step zero solution 

For the first moment of the test, the following output will be produced: 

(1) Basic wake information: 

(2) Solution iteration history: 

(3) Aerodynamic data for patch one which consists of a series of tables 
that lists for each of the panels in one chordwise row the x, y and z coordinate of the panel control 
point, as well as Cp, Vx, Vy, Vz, V, DUB, DPHI/DT, and MACH. This process is repeated for each 
row of panels of a patch. 

(4) Section parameters for the series of rows of panels which 
aerodynamic data were just listed. The section parameters are CL, CD, CY, Cm, Cn, and Cl. These 
values can be regarded as the two-dimensional values for that strip of panels. They are given in two 
reference sets: wind, and body axes. 

(5) Aerodynamic data and section parameters for all other patches . The 
same format as presented for the first patch is followed for all remaining patches. 

(6) Wind axis force and moment coefficients . In this section, total CL, 
CD, CY, and Cm are given for each of the patches, each component, and finally for the whole 
assembly. 

(7) Body axis force and moment coefficients . In this section, total CN, 
CA, CY, and Cm are given for each of the patches, each component, and finally for the whole 
assembly. 


c. Time-step 1 and above 

For the remaining time step except the last one, DATA6 will produce a condensed 
output in that it will not generate all the tables of aerodynamic data as for time step zero. Thus, all 
the steps listed above are executed except (3) and (5). 


d. The last time-step 

For the last time-step, all of the seven steps listed in time-step zero are executed 
again. The results will be different than when the flow just started at time step zero. 
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c. 


CONCLUSION 


The theoretical basis for PMARC, a low order potential-flow panel code for modelling 
complex three-dimensional geometries has been outlined. The input file was also detailed with the 
emphasis on the fact that although not all the fields apply to all cases, they all have to be filled with 
values for the software to run successfully. Finally, the output file was reviewed. Its format was 
studied along with the type of results that it holds. 
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III. SOFTWARE VALIDATION - CONVENTIONAL WINGS 


A. INTRODUCTION 

The objective of this chapter is to show that PMARC is a panel code that closely recreates 
most of the results from experimental testing, higher order panel codes, and theoretical analytical 
solutions. As it will be seen, there are restrictions to its capabilities in the area of velocity 
distribution over delta wings at zero angle of attack. This PMARC review will be done by modelling 
various types of wings for which analyses have been conducted and results are available. 

Specifically, two types of swept back wings, four delta wings, various rectangular wings, 
an F5 wing and a parabolic arc delta wing will be studied. The lift-curve slopes, the pressure 
coefScient distribution in both chordwise and spanwise direction, the velocity distribution and the 
equivalence rule results will all be comparison factors for our validation process. The following 
chapter will conduct a similar study on three bodies of revolutions. 

As mentioned in the previous chapter, PMARC operating instructions will also be presented 
here. Therefore the structure of the following sections, each corresponding to a wing configuration, 
will be as follows. A description of the wing studied will first be given. Then, an overview of the 
theoiy and results serving as the reference data will be presented, followed by the results obtained 
with PMARC, shown compared to the reference data, and complemented by explanatory comments 
and discussion. Also, the specific PMARC input steps for modelling the subject wing will be 
explained. Finally, the output data retrieval technique is presented. 

B. SWEPT BACK WING I 

A wing with relatively simple geometry was considered for the first comparison. The wing 
has an aspect ratio (AR) of 5, a taper ratio of unity, and an uncambered flat plate section. Since the 
taper ratio is unity, the leading edge, the quarter-chord line, the three-quarter-chord line, and the 
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trailing edge all have the same sweep, 45 degrees. All the dimensions and coordinates being in terms 
of the span b, an arbitrary value of b will be assumed at 10. The finite geometry of the wing is as 
shown in Figure 3.1. 



Figure 3.1 Four-Panel Representation of a Swept Planar Wing, Taper Ratio of 
Unity, AR=5, Sweep=45 Degrees. (From Ref. [4]) 
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1. Theory, Bacl^round and Reference Data 

Several lifting flow field computation methods have been developed. Probably one of the 
simplest is the Vortex Lattice Method (VLM). The VLM represents the wing as a planar surface on 
which a grid of vortex panels is superimposed. The velocities induced by each horseshoe vortex at 
a specified control point are calculated using the Bio-Savart law. Equation (3.1). 

dV = - 

47tr^ 

A summation is performed for all control points on the wing to produce a set of linear algebraic 
equations for the horseshoe vortex strengths that satisfy the boundary condition of no flow through 
the wing. The vortex strengths are related to the wing circulation and the pressure differential 
between the upper and lower wing surfaces. Berlin & Smith (Ref. [4]) present the derivation of the 
equations from the Bio-Savart law, and apply them to the case of the swept back flat plate shown 
above. The experimental results against which the PMARC and VLM results will be compared are 
reported in Ref. [5]. 

2. Results and Discussion 

As shown on Figure 3.2, the coefficient results obtained by PMARC correspond very closely 
to the ones calculated with the Vortex Lattice Method. However, although good for the low angles 
of attack, both sets seem to gradually increase their spread with the experimental results. That does 
not discredit PMARC or VLM, as this behavior was to be expected. PMARC and VLM are both 
linear methods that do not account for boundary layer separation and stall. That is why the slope of 
the experimental results shown on Figure 3.2 start to decrease as the stall angle approaches, while 
the other two methods keep showing a steady increase of lift coefficient with increase in angle of 
attack. 


r) 


(3.1) 


17 





Figure 3.2 Comparison of The Lift Coefficient as a Fvinction of Angle of Attack. 
Swept Planar Wing, Taper Ratio of Unity, AR=5, Sweep=45 Degrees. 


3. PMARC Input Instructions 

The following details complement Chapter n which presented the theory behind PMARC 
and its major operating instructions. This subsection will reiterate certain important details and 
introduce new factors that were intentionally left for this section for simplicity reasons. The 
following headings refer to the PMARC input file, DATA5, shown in Figure 3.3. 
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a. Basic Input Section 

When referring to the PMARC manual (Ref. [1]) for the definition of each of the 
variables, it becomes rather straightforward which values are to be entered in each field of the Basic 
Input Section. As mentioned in Chapter n, for most of the cases, only the following fields require 
attention: LENRUN, NTSTPS, DTSTEP, ALDEG, and SREF, although we will also continually 
fill CBAR, and SSPAN. Some others like VINF and VSOUND can be entered with actual values 
to represent the fi'eestream velocity and the altitude but these won't have any impact on the lift 
coefficients and pressure distribution results. It is to be noted that the lift, normal force, moment 
and drag coefficients are sensitive to the fields CBAR, SREF, and SSPAN, as the velocity 
components, the pressure coefficient, DUB, and dPHI/dT are not. The remaining fields should be 
filled with zeros or numbers taken fi'om existing DATA5 files. Although they too don't have any 
impact on the above mentioned coefficients, numbers have to be entered for the software to operate. 

One important point is the steady state status. Because of the time-stepping 
capability of PMARC, results are generated right from the first moment of the 'flight', while the 
circulation is initiated and the starting vortex generated. If the user is interested in the steady state 
results, which is when the starting vortex is far enough from the trailing edge (five or six chord 
length), a proper combination of number of time-steps (NTSTPS) and distance between each time- 
step (DTSTEP) will be required. Using the freestream velocity, this calculation is easy to perform. 
For instance, a wing with a chord of 11 units would require the starting vortex to be at a minimum 
of 55 units behind the trailing edge in order to assume steady state. If the freestream velocity is 11 
units/sec, then 5 seconds of flight are needed (55 units divided by 11 units/sec). Then any 
combination of NTSTPS and DTSTEP whose product is at least 5 units will ensure steady state. 

b. Assembly and Component Coordinate System Information 

All the fields but two can be copied from those shown in Figure 3.3 for most of the 
cases we will be concerned with. The two fields that may require special attention are NODEA and 
NODEC. Their values indicate respectively if another assembly coordinate or component coordinate 
will be defined. Such additional coordinate systems could be required when modelling a wing-body 
combination, or two wings alone. 
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c. Patch and Section Coordinate System Information 
All the fields in PATCH 1 are straight forward when referring to the PMARC manual 
(Ref. [1]). The point to note here is the existence of more than one PATCH 1 section within the input 
file. One introduces the patch that will be used to create the wing according to the geometric 
specifications given in the Section Coordinate section (SECTl and SECT2), and the second one is 
for the wing tip which is further defined in PATCH2. 

In SECTl, two fields are not relevant in steady state: ALF and THETA. All the 
others are crucial for a correct model definition. In this case, the origin of the root airfoil section 
definition is at coordinates (0,0,0). The root chord is two units. INMODE was set at 5 because in 
SECT2 we will generate a NACA 4 digit airfoil section for which PMARC has automatic section 
coordinate generation routines. All that is required are the thickness to chord ratio (RTC), the 
maximum camber to chord ratio (RMC), and the chordwise position of the airfoil's maximum 
camber (RPC). In this case, we are dealing with a flat plate. These three combinations all produce 
the same results: 

RTC=0.01 RTC=0.01 RTC=0.01 

RMC=0 RMC=0 RMC=0 

RPC=0 RPC=1 RPC=0.5 

In our case, the second approach was taken which corresponds to the NACA designation of 0001. 
It becomes obvious why the three combinations shown above work by noting that for a flat plate, 
there is no camber (RMC), therefore the position of maximum camber as described by the input field 
RPC has no impact, no matter what value is entered. 

The second SECTl section defines the location and section of the wing tip, or more 
specifically that flat surface where the wing tip would be attached. We see that in our case, the 
origin of that last patch's section is at the coordinates (5,5,0), and that its chord is also 2. As defined 
in Ref. [1], INMODE and TNODS seen in SECTl and PATCH2 respectively are important fields. 
INMODE set to zero in the second SECTl tells PMARC to copy the section definition of the 
previous SECTl. The TNODS fields indicate the status of the patch being defined: first, 
intermediate, last, etc. 
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Ten panels per side (TNPS=10) are requested spanwise for the specified patch. That 
means 20 panels spanwise for the complete wing on each surface. Fifteen panels are prescribed 
chordwise (TNPC=15) for the specified patch. That transpose into 15 panels on each surface of the 
wing. Note that TNPS can be specified only in the second SECTl section. If both TNPS were set 
to zero, PMARC would use the points of the section definition it created to determine the number 
of panels. 


As mentioned earlier, the second PATCH 1 section along with PATCH2 are for the 
wing tip, and a circular arc tip patch was requested (ITYP=2). 

Two fields were not discussed yet but become very important for the accuracy of 
the results. They are TINTS and TINTC. They are found in SECTl, SECT2, and PATCH2. They 
define the type of panel spacing on each patch (or wing and body surfaces). As it will be seen later 
for other configurations, the user might want to distribute more panels in areas where a sharper 
gradient of a parameter is known to exist. Figure 3.4 shows the available options of panel spacing 
for an arbitrary wing. 


d. Wake and Streamlines 

The streamlines are not required for our studies. However, the wake should be 
turned on (IDWAK=1) in all wing cases as it affects the flow field all around the wing. For the 
simple cases studied in this thesis, the type of wake selected (rigid or flexible) in the IFLXW field 
is not a crucial decision as this will not change the results. 

4. PMARC Output - Data Retrieval 

The results we are interested in are the lift coefficients (Cj,) for various angles of attack. In 
order to do so, our input file has to be run repetitively, each time modifying the angle of attack. This 
will create a series of output files, DATA6, from which the will be read. Several lift coefficients 
are presented in this output file: assembly, component, and total, in term of wind or body axes. For 
the present case, we need the total wing lift coefficient in terms of the wind axis which refers to the 
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freestream. This value is found towards the veiy end of the file under the heading Wind Axes - Total 
coefficients. The lift coefficient corresponding to our case is labelled CL. Note that when 
streamlines are requested, the related results are placed at the end of the file, following the 
coefficient results. 



Figure 3.4 PMARC Options of Chordwise and Spanwise Panel 
Distributions. 











c. 


SWEPT BACK WING II 


Our second project was to model a cambered swept back wing. The wing considered had 
a NACA 64A010 section, an AR of 3, a tip to root chord ratio of 0.5 and a 45 degrees swept back 
of the quarter-chord line. A value of the root chord Cr was arbitrarily set at 10, allowing the finite 
coordinates as shown on Figure 3.5. 



Figure 3.5 Swept Back Wing, NACA 64A010 Section, AR=3, Taper 
Ratio=0.5, Sweep=45 Degrees. 


1. Theory, Background, and Reference Data 

Margason et al. (Ref. [6]) compared computed aerodynamic coefficients using the Vortex 
Lattice Method, source panel method, low-order surface potential distributions, and high-order 
surface potential distributions to experimental data. This provided an ideal opportunity to see how 
PMARC would measure up. 
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2 . 


Results and Discussion 


The results are shown on Figure 3.6. It can be observed that all surface potential and source 
panel methods overpredict the experimental data. However, each of them show little difference 
between their lift coefficients, PMARC included. The results for the Vortex Lattice Method are not 
presented, but Margason reports that they predict the experimental data very well, due to the fact that 
VLM neglects both thickness and viscosity effects. 



Figure 3.6 Comparison of The Lift Coefficient as a Function of Angle of Attack. 
Aft Swept Wing, NACA 64AO 10 Section, AR=3, Taper Ratio=0.5, Sweep=45 
Degrees. 


Although PMARC did not recreate exactly the experimental data for the NACA 64AO 10 
swept back wing, its deviation from the latter is still acceptable. Furthermore, its results compare 
very well with other computer methods, even with panel codes of higher order. 
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3. 


PMARC Input Instructions 


The following sub-subsections will not repeat information presented for the swept back flat 
plate that still hold for the cambered swept-back wing. Instead, we will concentrate on the 
differences, the particularities, and some lessons learned. The input file is shown in Figure 3.7. 

a. Basic Input Section 

All the fields in this section can remain the same as for the flat plate case, except 
for CBAR, SREF, and SSPAN. Of course, the angle of attack (ALDEG) has to be adjusted as 
required. Also, VINE and VSOUND could be changed, although there is no need to do so. 

b. Assembly and Component Coordinate System Information 

No changes are required in these two sections. As it can be seen in the definition 
section of Ref. [1], these terms become important when multiple components are involved, or when 
a component needs to be rotated or translated. 

c. Patch and Section Coordinate System Information 

The PATCHl sections for both the flat plate and cambered cases are identical. The 
PATCH2 lines have for only difference the panel field (TNPS). It could have been the same as in 
each case it was an arbitrary decision. This field defines the number of chordwise (longitudinal) cuts 
on the wing tip, which along with the vertical cuts will define the number of panels on the tip (Figure 
3.8). TNPS of PATCH2 can be kept low without altering the results. Fewer panels speeds up the 
computation and in this case alleviates a message in the output file warning the user that some panels 
have virtually no area. Although this warning did not seem to affect any of the results, caution 
should be exercised and more investigation should be pursued in this area. It is to be noted that 
Figure 3.8 shows a flat tip only for panels illustration purposes. All our wings have rounded tips 
(except in the case of the delta wings). 


26 



Q Q Q Q Q _ as 
z as z z z o u 
u u u u u z 


I o ^ • ~:c o o o • o 
iiooo^ • *.oi 
>x«i •a.ouooKfl 
S 2 oxewi zi a x< 
HiiaaMNOiujSz 
SMMaCK &NUDS 07 


I --^oaoeoes - »o «. ^ 
bOIttX ‘^OWIOO. 

k]lop2oxoi Runi lu 

ao.=><w 5 -x>-wojgo 2 yw 
r<OOXXtf>it* 2 SSeaO 0 
:) DC u> H X £3 » <0 2 z u z z 


vaina iooo*«krn 

cucaw-i-^ax •■ •i^«^fr<szooi 
ZWH( 4 « lOaoXOR IMOOa »H 
)-iVM|4Xk>MZ a < * 2xv3a<sk,aL.ZM 

e<Sxw>zQi-tx 2 xao,oSN 2 yaJ 

ij»a 2 za:>ao.SQXu 2 zzzzxz 


I u o a 

>HZiaoeoooooc 

U<<H XlbXOOOOOOC 
£o.X[-<>-3Zooooooc 
MW m << M. 


Figure 3.7 PMARC Input File for an Aft Swept Wing, NACA 64AO 10 Section, AR=3, 
Taper Ratio=0.5, Sweep=45 Degrees 
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Figure 3.7 Cont'd PMARC Input File for an Aft Swept Wing, NACA 64A010 Section, 
AR=3, Taper Ratio=0.5, Sweep=45 Degrees 
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(From Ref. [1]) 

The SECTl and SECT2 combination is different for this case of cambered wing. 
Instead of a four digit NACA airfoil section, we now have a section definition for which PMARC 
has no automatic construction routine. We then have to enter the coordinates of the airfoil section 
definition. For this reason, SECT2 is not required anymore, as it is only used in the case where a four 
digit NACA airfoil is wanted. 

We will still observe two SECTl sections: one for each end of the patch (or wing). 
The differences between the two reside in the coordinates of the origins of each section starting and 
ending the patch. Also, the scaling is different, meaning that the outside edge of the patch (wing tip) 
has a different length than the inside edge (wing root). An important field here is INMODE. It is set 
at four in the first SECTl of the cambered wing input file to indicate that the patch coordinates 
follow. In the flat plate input file, the first INMODE was set at five requesting that PMARC generate 
an NACA four digit airfoil section. The second INMODE is set at zero, instructing PMARC to copy 
the section definition of the previous section. The same thing could have been done with the 
cambered wing instead of redefining the section definition for the tip. 
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The manual entry of the section definition should be done carefully. Some entry will 
not cause a system error, but will produce bad results. For instance, if the patch coordinates are 
entered starting at the origin (0 0 0), PMARC will try to satisfy the Kutta condition at the leading 
edge instead of at the trailing edge. Therefore, an important lesson learned while exploring PMARC 
was to understand how the patches are defined and folded (Ref. [1]). For the case of any wing, the 
section definition always starts at the trailing edge, going under the wing, around the leading edge 
and up towards the trailing edge (Figure 3.9). 



(From Ref. [1]) 
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There is more than one way to use the coordinate system to define a section. That 
is controlled by the field INMODE in SECTl. Depending on its value, the three columns of data will 
have different meanings (refer to page 53 of Ref [1]). In the present case, the three columns 
respectively represent the x, y, and z coordinates of each point defining the section. This section 
definition of the NACA 64A010 was taken in Ref. [7]. 

As well illustrated in Ref. [2], and in the test case of Ref. [1], more than one means 
are available to shape the required dimensions into our wing sections. In the present case, we elected 
to non-dimensionlize the section definition and then scale it to the chord and tip dimensions 
respectively. Another way would have been to enter the global coordinates of the root and tip 
sections. 


d. Wake and Streamlines 

As mentioned in the previous section, we need to leave the wake active (IDWAK=1) 
as we do not need the streamlines option. 

4. PMARC Output - Data Retrieval 

Once again, as for the swept back wing of the previous section, we are interested in are the 
lift coefficients (Q) for various angles of attack. In order to do so, our input file has to be run 
repetitively, each time modifying the angle of attack. This will create a series of output files, 
DATA6, from which the Q will be read. Several lift coefficients are presented in this output file: 
assembly, component, and total, in term of wind or body axes. For the present case, we need the 
total wing lift coefficient in terms of the wind axis which refers to the freestream. This value is 
found towards the very end of the file under the heading Wind Axes - Total coefficients. The lift 
coefficient corresponding to our case is labelled CL. Note that when streamlines are requested, the 
related results are placed at the end of the file, following the coefficient results. 


31 




D. DELTA WINGS 


A good amount of experimental data is available for delta wings, allowing a wider spectrum 
of comparison. Four different delta wings were modeled, each having different aspect ratio (AR), 
but all with aNACA 0012 airfoil section designation. In order to obtain finite geometric parameters, 
the span was arbitrarily set at 10 in all cases. Figure 3.10 shows those configurations. 



1. Theory, Background and Reference Data 

Prediction and control of flow separation is a major aerodynamic consideration in wing 
design. Effort to avoid undesirable compressibility effects leads to greater sweep angle and reduced 
section thickness. However, these efforts are offset by the increasing difficulty to prevent boundaiy- 


32 




layer separation. Schlichting and Truckenbrodt (Ref. [8]), and also Ashley and Landahl (Ref [9]) 
studied the influence of the delta wing aspect ratio on the lift coefficients and the boundary layer 
separation. The purpose of our study is not to draw any conclusion with regard to flow separation 
and AR influence using our PMARC computation but rather to try to recreate some of the 
experimental results obtained by the above mentioned researchers. 

Our comparative study will present five sets of data. We will first compare the lift 
coefficients of the delta wings of AR 0.83 and 3.97 obtained by PMARC with the available 
experimental data. Then a plot of the lift curve slope as a function of the AR will be built with the 
PMARC results of four different AR wings. Finally, the chordwise and spanwise pressure 
distributions for the delta wing of AR=3 will be established using PMARC and compared to the 
literature. 


2. Results and Discussion 

The first delta wing modeled had an aspect ratio of 0.83. The lift coefficient for various 
angles of attack as obtained with PMARC is shown compared to the experimental data on Figure 
3.11. Very good agreement exists at low angles of attack but above fO degrees, PMARC 
underpredicted the experimental data. This situation was encountered before in this paper and as it 
was mentioned then, this disparity at higher angles of attack does not invalidate PMARC. These 
results were expected as PMARC only models linear angle of attack effects and therefore does not 
account for nonlinear phenomena encountered at high angles of attack. Therefore, for the lower 
angles of attack, PMARC recreates the experimental data very well. 

The same phenomena occur with the delta wing of aspect ratio 3.97 where the CL obtained 
by PMARC steadily and linearly increases with angle of attack (Figure 3.12). The difference here 
is that PMARC overpredicted the lift coefficient values. Again, for the linear portion of the 
experimental data, PMARC is in close agreement. 
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Angle of attack (Degrees) 

Figure 3.11 Comparison of The Lift Coefficient as a Function of Angle of 
Attack. Delta Wing, NACA 0012 Section, AR=0.83. 



Figure 3.12 Comparison of The Lift Coefficient as a Function of Angle of 
Attack. Delta Wing, NACA 0012 Section, AR=3.97. 
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Lift-Curve Slope, dCI/d(atpha) - /radian 


The next figure (Figure 3.13) shows the lift curve slope as a ftinction of aspect ratio. Again, 
the PMARC results display good agreement with past work (Ref. [9]). It is shown there that the 
lifting-surface theory predicted very well two sets of experimental data produced at different time 
by different people. We superposed the PMARC output onto these combined results to show the 
very good agreement. 

We have so far observed that PMARC seems to compute fairly accurately the lift 
coefficient. Another parameter should now be introduced in order to better assess this panel code. 
As such, the last two sets of results presented in this section will deal with the pressure distribution 
for a delta wing of aspect ratio equal to three and an angle of attack of 10 degrees. 



Figure 3.13 Comparison of Lift-Curve Slope as a Function of Aspect Ratio. 
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The chordwise and spanwise distributions are shown on Figures 3.14 and 3.15. Figures 
3.16 and 3.17 taken from Ref. [17] show experimental spanwise pressure distributions. Some 
irregularities appear to be present in the PMARC output of the spanwise distribution towards the 
edges of the spanwise sections considered. The pressure do not seem to come back down as shown 
on Figure 3.17. Either we do not have enough points or the panel layout can not pick up the correct 
behavior of the flow in these regions of greater gradient. Note that the PMARC spanwise 
distribution plotted are in reality approximations to the real spanwise distribution. What is plotted 
is actually the pressure distribution along a slanted spanwise cross-section that corresponds to panel 
munbers easy to retrieve out of the output file. More detail on this topic will be given in subsection 
4. Also, the reader should keep in mind that Figures 3.16 and 3.17 represent visous flows. 



Figure 3.14 Chordwise Pressure Distribution at 10 Degrees Angle of Attack. Delta 
Wing, NACA 0012 Section, AR=3. 
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Figure 3.15 Approximate Top Spanwise Pressure Distribution at 10 Degrees Angle 
of Attack for 1/3 and 2/3 Chord Locations. Delta Wing, NACA 0012 Section, AR=3. 



Figure 3.16 Spanwise Pressure Distribution Over a Flat and Cambered Delta 
Wing, AR=1.15, AOA=25 Degrees. (From Ref [17]) 
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3. PMARC Input Instructions 

Five input files were required to obtain all the results presented in this section. They were 
all identical except for the section sizing and coordinates (SECTl) and SREF. These sets of fields 
correspond to the specific wing geometiy as well as to the root and tip section origins. One sample 
for the delta wing of aspect ratio 0.83 is shown in Figure 3.18. The fields mentioned above are also 
the only differences between the delta wing input files and the swept back flat plate input file of our 
first test case. Therefore, no time will be spent here to go over each section of the input file. Instead, 
we will discuss a particularity of the delta wing modelling in PMARC. Also a sensitivity test 
performed on some fields will be discussed. 
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a. Shaping a delta wing tip 

In order to bring our NACA 0012 airfoil section from its full chord length at the root 
to a point at the tip, an initial attempt was made by scaling it to zero (SCALE=0); but the system 
crashed. To alleviate this situation, it was found that a number other than zero had to be entered. So 
to keep the integrity of oiu' geometry a very small number was selected, SCALE=0.05. In fact, that 
was the smallest value that would allow the system to run successfully. This number could vary as 
it will be seen in a later case where it was successfully set at 0.01. This variation is a function of the 
overall geometry of the wing, and may be related to the number of panels on the surface. 

Although it has not been confirmed yet, the above situation may be related to the 
size of the panel in the region of the tip. Panels that are too small may cause problems. This area 
needs more investigation as we did run cases where the output warned the user about panels having 
virtually no area, but still produced valid lift coefficients and pressure distributions. 

b. Testing Certain Fields 

The software sensitivity to the number of panels on the tip was brought up before, 
but we want to warn potential PMARC users that it is even more true in the case of very small tips 
like for the delta wing. The wing tip section of the input file, the TNPS field of the PATCH2 line 
should be set at a very low value. The smaller the tip, the smaller the value. In the case of the delta 
wing, TNPS=1 seemed to be the most efficient value for a successful run. 

4. PMARC Output - Data Retrieval 

We already discussed where to find the lift coefficient in the output file DATA6. We need 
now to indicate how to obtain the data for the pressure distribution shown in the above plots. It will 
be different depending along which direction (spanwise or chordwise) we need the distribution. 
Either way, a clear picture of how the panels are labelled is crucial. Figure 3.19 details this 
important information for an arbitrary case where TNPS is 8 and TNPC is 12. 
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Figure 3.19 PMARC Panel Labelling Around a Delta Wing. 


41 






a. Chordwise Pressure Distribution 

In order to plot the steady state chordwise pressure distribution, the user needs to 
scroll down the output file until the time-step corresponding to steady state is reached. If the time- 
step fields NTSTPS and DTSTEP are filled in accordance with the instructions provided in Chapter 
n, the data corresponding to the last time-step (towards the bottom of the file) are within the steady 
state regime. In the present case we wanted to plot the distribution at the root, and the output labels 
the root panels fi-om one to thirty under column one. 

As will be seen later, chordwise pressure distribution at spanwise locations other 
than the root can also be plotted using the PMARC output. The only additional step is to locate the 
column of data with a constant y value corresponding to the spanwise location desired. 

b. Spanwise Pressure Distribution 

In order to precisely retrieve the Cp value for each spanwise station along a 
spanwise cross-section, one has to select one panel per column whose coordinates are the closest to 
the cross-section and station under consideration. This could be a lengthy process, but some 
computer program can be written to do it automatically. What was done with the delta wing 
spanwise pressure distribution was an approximation believed to be good enough to analyze the 
trend. 


The fifth to last panel from each of the ten columns was selected and the Cp value 
recorded for each. The process was repeated with the tenth from last panel of each column. As the 
panel selected in each column are closer to the trailing edge, the panels form a line more and more 
parallel with a normal spanwise distribution. As the panels get away from the trailing edge, this line 
moves more and more towards the leading edge, ultimately merging with it. Consequently, the 
results shown for the 2/3 chord location are more accurate than those of the 1/3 chord location, but 
the trend is reasonably the same. 
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E. RECTANGULAR WINGS 


The rectangular wing is the simplest from the point of view of structural design and analysis. 
In this section we will put PMARC to the test with rectangular wings of aspect ratio 1 , 4.5, and 2, 
respectively. All wings will have a NACA 0012 cross-section, except for a study shown at the end 
where the wing of AR=7 is gradually thinned out to observe the impact on the pressure distribution. 
Figure 3.20 shows the relative size of the three wings. For the PMARC input file, we arbitrarily 
selected all chord lengths to be 1, leading to a span of 7,4.5, and 2 respectively. 
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Figure 3.20 Rectangular Wings of Various Aspect Ratio. 


1. Theory, Background and Reference Data 

Two main references were consulted in order to compare the results obtained with PMARC. 
The first one is a Vortex-Lattice software output produced in 1971 by Margason and Lamar (Ref 
[10]). The program represents the lifting planform with a vortex lattice and solves the system of 
equations in a manner summarized in section ni-B-l and detailed in the above mentioned reference. 
The Vortex Lattice Method is not the most accurate method for estimating subsonic aerod 5 mamic 
characteristics of lifting planforms. Consequently, the validation of the PMARC results should 
account for this fact. 
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The second reference is a two-dimensional computer panel code named PANEL developed 
by Dr. M. Platzer (Ref. [11]). For any airfoil coordinates and a selected number of panels, it 
computes the major aerodynamic coefficients. We will use PANEL pressure distribution results for 
a NACA 0012 to compare with the PMARC output of a large aspect ratio (20) NACA 0012 
rectangular wing. Although PMARC is a three-dimensional code, the mid-span results for a high- 
aspect-ratio wing can be regarded as the same as the two-dimensional values as the wing tips have 
negligible effect at the root. 

2. Results and Discussion 

a. Comparing with the Vortex-Lattice Results 

In addition to comparing the lift-curve slope of PMARC and Ref [10], the following 
results will show the influence of the number of panels on the final output. Figure 3.21 presents the 
results obtained by Margason & Lamar superposed by those computed by PMARC. For all three 



Figure 3.21 Effect of Panel Arrangement on Lift-Curve Slope for Rectangular Wings 
of Various Apsect Ratio. 
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wings, the lift-curve slopes as obtained by both methods asymptotically decrease as the panel density 
increases. In all cases PMARC overpredicted the Vortex-Lattice results by five to 10 percent, the 
largest disparities belonging to the smallest aspect ratio. Also shown on the figure is the result when 
the chordwise panel density is significantly higher than what is presented. It is suspected that the 
latter results are the most accurate simply because more panels allow for better wing definition. But 
again, these results do not differ from Margason and Lamar's highest panel density results by more 
than 7.7%. 


b. Comparing Chordwise Cp Distribution with the 2-D Panel Code 

In order to efficiently use PMARC to get exact two-dimensional results, a wing of 
infinite span should be constructed. As discussed earlier, a good approximation is a wing of 
sufficiently large aspect ratio. For the case of the rectangular wing, we elected to test a 20 unit wing 
span with an aspect ratio of 20. The object of the tests was to vary the thickness of the airfoil section 
from that of a NACA 0012 to a NACA 0001, and to compare the pressure distribution of each wing 
at an angle of attack of ten degrees, as computed by PMARC and PANEL. 

As shown on Figures 3.22 to 3.28, very good agreement between the two codes 
exists for all airfoils. It is to be noted that the lack of smoothness in the PMARC Cp curves is due 
to the small number of panels used in the computation. That also explains why the peaks are 
sometimes cut-off. Additional panels would have meant more data points resulting in a better 
resolution of the curve. 

Figure 3.29 is a composite Cp plot for all seven wings as computed by PMARC. 
This plot was included for general information to show on a relative scale how the pressure 
distribution changes as the wing gets thinner. 
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Figure 3.22 Comparison of The Pressure Coefficient as a Function of Chord 
Location. Rectangular Wing, AR=20, AOA=10 Degrees, NACA 0012 Section. 



Figure 3.23 Comparison of The Pressure Coefficient as a Function of Chord 
Location. Rectangular Wing, AR=20, AOA=10 Degrees, NACA 0010 Section. 
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Figure 3.24 Comparison of The Pressure Coefficient as a Function of Chord 
Location. Rectangular Wing, AR=20, AOA=10 Degrees, NACA 0008 Section. 



Figure 3.25 Comparison of The Pressure Coefficient as a Function of Chord 
Location. Rectangular Wing, AR=20, AOA=10 Degrees, NACA 0006 Section. 
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Figure 3.26 Comparison of The Pressure Coefficient as a Function of Chord Location. 
Rectangular Wing, AR=20, AOA=10 Degrees, NACA 0004 Section. 



Figure 3.27 Comparison of The Pressure Coefficient as a Function of Chord 
Location. Rectangular Wing, AR=20, AOA=10 Degrees, NACA 0002 Section. 
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Figure 3.28 Comparison of The Pressure Coefficient as a Function of Chord Location. 
Rectangular Wing, AR=20, AOA=10 Degrees, NACA 0001 Section. 
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Figure 3.29 PMARC Chordwise Pressure Distribution for Various NACA 
Rectangular Wing Sections, AR=20, AOA=10 Degrees. 
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c. Spanwise Cp Distribution 

Two wing thicknesses were tested in this section for the rectangular planform of 
aspect ratio 20. Figure 3.30 shows the comparative pressure distributions for the NACA 0001 and 
0012, taken at one third of the chord length of each wing. Figures 3.31 and 3.32 show each wing 
spanwise Cp individually but taken at two different chordwise locations. 



Figure 3.30 Spanwise Top Pressure Distribution at 1/3 Chord for Two NACA Sections. 
Rectangular Wing, AR=20, AOA=10 Degrees. 


Although we do not present experimental data to corroborate the PMARC spanwise 
pressure distribution results, an intuitive analysis leads to a positively conclusive validation of the 
PMARC output. First, in all cases PMARC produced the expected elliptical distribution. Then, in 
Figure 3.30, the plots show higher Cp's for thicker wings, which was also to be expected. Figures 
3.31 and 3.32 again present the expected behavior. However, one area of concern remains towards 
the tips where a sudden peak indicates that the panel layout may not pick up the correct behavior in 
these area of greater gradient. 
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Figure 3.31 Spanwise Top Pressure Distribution at 1/3 and 2/3 Chord Locations. 
Rectangular Wing, AR=20, AOA=10 Degrees, NACA 0012 Section. 
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Figure 3.32 Spanwise Top Pressure Distribution at 1/3 and 2/3 Chord Locations. 
Rectangular Wing, AR=20, AOA=10 Degrees, NACA 0001 Section. 
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3. 


PMARC Input Instructions 


Except for the fields specific to the wing geometry (CBAR, SREF, SSPAN, SCALE, SIX, 
and STY), the test and panel prescription (ALDEG, TNPS, and TMPC) nothing is changed from the 
previous test case input file, the delta wing, to the rectangular wings of various aspect ratios. Figure 
3.34 shows the input file for the rectangular wing of aspect ratio 7. Similarly, the input files of the 
other aspect ratios studied vary only by the above fields. 

4. PMARC Output - Data Retrieval 

As mentioned in the previous section, the data retrieval required in the PMARC output file 
when interested in plotting pressure distributions necessitates a complete understanding of the panel 
numbering method. Figure 3.33 illustrates this protocol for the case of the rectangular wing. The 
data retrieval for the lift coefficient, the chordwise and spanwise pressure distribution was discussed 
in the previous section. 



Figure 3.33 PMARC Panel Layout for a Rectangular Wing and Its Tips. 
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F. F5 WING 

The F5 wing is the subject of an independent study. The results of the latter will be 
combined and analyzed with those of this paper in subsequent work. Also, some unsteady 
aerodynamic results will be introduced later in this paper. Here we merely present the PMARC 
predictions for steady flow. Figure 3.36 shows the half wing non-dimensionalized coordinates. 
From these coordinates, all the geometry was defined, an aspect ratio of 2.98 calculated, and a 
thickness ratio of 0.0480 derived. Figure 3.35 is the PMARC version of the F5 wing. The airfoil 
cross-section corresponds to a NACA 65004.8 whose coordinates were provided by Ref. [13]. The 
leading and trailing edge points were not presented in that reference and therefore had to be 
introduced using a curve fit of the available points. Figure 3.37 shows the given and the fitted points 
of the wing section. 

The dimensionless chord coordinates exceed unity on Figure 3.37 to show the relative 
magnitude of the appended edges. However all data points were divided by the new chord length 
in the PMARC section definition in order to conserve a chord length of unity. 



Figure 3.35 PMARC Representation of The F5 Wing. 
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Figure 3.37 Leading and Trailing Edge Fitting. F5 Wing Cross-Section. 


1. Theory, Background and Reference Data 

a. Pressure Distribution - Experimental Results 

The data used for comparison with the PMARC results are from experiments 
presented in Ref. [12]. Figure 3.36 displays the wing layout for the experiment with the pressure 
orifices, in-situ transducers, and accelerometers (Kulite and Endevco). One of the results presented 
in this study is the pressure distribution at different spanwise locations at zero degree geometric 
angle of attack. The experimental results considered were taken at a Mach speed of 0.6 and 
accoimted for compressibility effects. As will be shown, the Prandlt-Glauert compressibility 
correction factor will be applied to the PMARC incompressible results in order to account for 
compressibility effects. 

b. Velocity Distribution - Area Rule Analytical Results 

In the early days of supersonic flights, aircraft designed to penetrate the sound 
barrier ran into difficulty as the thrust of their engines could not overcome the large peak drag near 
Mach 1. It is after the advent of the area rule that dramatic breakthrough in high-speed flight started 
to flourish. Richard T. Whitcomb, the father of the area rule, reasoned that the variation of cross- 
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sectional area for an airplane should be smooth, with no discontinuities. That is how the famous 
Convair F-102 finally broke the sound barrier; after its fuselage was modified to allow a reduction 
in its cross-sectional area in the region where the wings were located, in such a way as to smoothen 
the total cross-sectional area distribution. 

The analytical derivations behind the area rule is presented later in Chapter HI, 
section G where we verified the area rule by comparison with Ref [15]. We will see there that the 
flow can be represented by a cross-sectional flow (or cross-flow) to which is added a spatial 
influence component. The velocity potential of the cross-sectional flow and the spatial influence are 
respectively expressed as follows: 


<P = 


lit 


lll[ 




_1_ 1 
2nJ-s(x)^ ^-z ^ ^ 


(3.2) 



QJx)[/n2x+ln2(x^-x)] + 


1 

471^0 |x-^| 


(3.3) 


Chapter lEL, section G presents various forms of the velocity components (velocity fields, or partial 
derivatives) of the above equations. 

In summary, the cross-sectional flow is produced in each y-z plane (Figure 3.38) 
by a distribution of two-dimensional sources along a spanwise axis. Each cross-flow is completely 
independent of conditions at any other y-z plane. 
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Figure 3.38 Low-Aspect-Ratio Wing and its Cross-Flow Potential. 

(From Ref. [14]) 



Figure 3.39 Substitution of Source Panels. (From Ref. [14]) 
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The spatial influence expression accounts for those parts of the wing situated farther 
upstream (and downstream in subsonic flow) from the y-z plane under consideration (Figure 3.39). 
The essential aspect of this theory is that the spatial term is a function only of the total cross- 
sectional area distribution with x. The shape of the wing (or wing-body combination) is irrelevant, 
only the cross-section area is important. This will become a vital characteristic when complicated 
shapes are to be modeled by computer software. An equivalent area body of revolution could be 
modelled instead which is the easiest geometry to build. Of course, the cross-flow still has to be 
solved. 


In the rest of this chapter, we will partially apply the concept of the two-component 
velocity potential on the F5 wing introduced earlier. We will discover that PMARC does not 
recreate valid potential distributions. For completeness, the results of the area rule using 
Equations (3.2) and (3.3) for the triangular wing will be presented in the next section. It will be 
shown that once more, PMARC can not reproduce an adequate velocity potential distribution. 

2. Results and Discussion 

a. Pressure Distribution 

Figure 3.40 presents the chordwise pressure distribution on the F5 wing when at an 
angle of attack of zero degree. The spanwise location that produced this distribution is very close 
to the root. It corresponds to the control points of the first row of panels whose edges define the 
root. Although not shown, these results agree with the experimental results. 

The next set of figures (Figures 3.41 to 3.44) better represents the agreement 
between PMARC and the experimental results. For four spanwise locations of the wing, the pressure 
distribution as calculated by PMARC is plotted against the experimental results. Also, since the 
latter account for compressibility and the former does not, the PMARC outputs were corrected with 
the Prandlt-Glauert factor and the results superposed on the respective figures. 
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Figure 3.40 Pressure Distribution as a Function of The Root Chord Location. F5 Wing, 
NACA 65A004.8 Section, AR=2.98, AOA=0 Degree. 


The PMARC output for the F5 wing did not allow for spanwise definition such as 
to have perfect correspondence with the spanwise locations selected for the experiment. However, 
the disparity is sufficiently small to be considered negligible in all but the first case where some care 
should be exercised as there is 11% difference between the PMARC and the experimental spanwise 
location. 


With this observation in mind, the reader can still appreciate that all the results 
display correspondent behavior, and for the last three figures especially, the Prandlt-Glauert 
correction seems to bring the PMARC results progressively closer to the experimental results as we 
approach the wing tip. 
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Figure 3.41 Comparison of The Chordwise Pressure Coefficient at a Spanwise Area. F5 
Wing, NACA 65A004.8, AR=2.98, AOA=0 Degree. 
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Figure 3.42 Comparison of The Chordwise Pressure Coefficient at a Spanwise Area. F5 
Wing, NACA 65A004.8, AR=2.98, AOA=0 Degree. 












Chordwise Location - Percent of Local Chord 

Figure 3.43 Comparison of The Chordwise Pressure Coefficient at a Spanwise 
Area. F5 Wing, NACA 65A004.8, AR=2.98, AOA=0 Degree. 



Chordwise Location - Percent of Local Chord 

Figure 3.44 Comparison of The Chordwise Pressure Coefficient at a Spanwise Area. F5 
Wing, NACA 65A004.8, AR=2.98, AOA=0 Degree. 
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Another result presented is the spanwise pressure distribution. Although we had no 
comparative data for these results, an intuitive analysis seems to validate them. Indeed, the general 
behavior combines those of the rectangular and delta wings. Spanwise cuts of the F5 wing in the 
region between the forward apex and the tip leading edge and similarly between the tip trailing edge 
and the aft apex correspond to cuts on a delta wing for which we see a bucket shape Cp distribution 
(Figure 3.45). Cuts in the region between the wing tips correspond to cuts on a rectangular wing for 
which we observe an elliptical Cp distribution (Figure 3.46). The PMARC spanwise pressure 
distribution on the F5 wing at zero degree angle of attack combines these behaviors and is shown 
on Figure 3.47. 



Figure 3.45 Typical Spanwise Pressure 
Distribution for a Delta Wing. 



Figure 3.46 Typical Spanwise Pressure Distribution 
for a Rectangular Wing. 
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Figure 3.47 Approximate Pressure Distribution as a Function of Spanwise Location. F5 
Wing, NACA 65A004.8, AR=2.98, AOA=0 Degree. 


b. Velocity Potential Distribution 

Comparison with the analytical results of Ref. [15] shows that the PMARC velocity 
potential output is erroneous. Figure 3.48 is the PMARC distribution of the x-component 
perturbation velocity u for various spanwise locations. It shows that the flow is faster as the 
spanwise location gets closer to the tip (that can be seen by noting that the peak velocities are higlier 
as z increases). Here, z is the spanwise location. (This notation does not correspond to the PMARC 
one but rather matches the technical note from which the area rule was taken (Ref. [15])). The 
perturbation velocity u is computed as shown in Equation (3.4). 


X °° 
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Figure 3.48 Perturbation Velocity on The Top Surface at Various Spanwise Locations. 
F5 Wing, NACA 65A004.8 Section, AR=2.98, AOA=0 Degree. 


This behavior is obviously not realistic. As the wing thickness tapers down toward 
the tips, the velocity perturbation should decrease towards zero. Figure 3.49 shows the distribution 
of the same parameters but for the lower side of the wing. Once again, the velocity peaks are 
wrongfully increasing toward the tip. 

Although we do not have any comparative data on the F5 wing spatial influence 
portion of the total flow, it was analytically solved using one form of Equation (3.3) shown later as 
Equation (3.23), and the results are plotted in Figure 3.50 for future consideration. The solved 
analytical steps can be found in Appendix A. 
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Figure 3.49 Perturbation Velocity on The Bottom Surface at Various Spanwise 
Locations. F5 Wing, NACA 65A004.8 Section, AR=2.98, AOA=0 Degree. 


Spatial Influonce 



Figure 3.50 Spatial Influence u-Component of The Velocity Potential. 
F5 Wing, NACA 65A004.8 Section, AR=2.98, AOA=0 Degree. 
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The cross-sectional portion of the total flow (Equation (3.2)) could not be solved 
because one of its terms, the two-dimensional cross-sectional segment area P(x,z), could not be 
constructed. Only unidimensional expressions could be derived with the numerical data available. 

3. PMARC Input Instructions 

At this point, the reader should have acquired enough knowledge about PMARC input to 
adapt it to various configurations. The F5 wing does not show any differences from what we have 
seen before, except for the usual basic test condition and geometry fields ( Figure 3.51). 

One aspect of consideration which has already been mentioned is the importance of the 
order in which the airfoil section coordinates are defined. Figure 3.52 shows a pressure distribution 
calculated by PMARC for the same F5 wing, but with a section definition starting at the leading 
edge instead of the trailing edge (Figure 3.53) as required by the PMARC protocol for a wing. As 
it can be seen, although some of the results are relatively close from the valid distribution, major 
irregularities exist, as the Kutta condition is not met at the trailing edge. (Note that on Figure 3.52 
the chord coordinates are not non-dimensionalized.) 

4. PMARC Output - Data Retrieval 

Again, all techniques for retrieval of the chordwise and spanwise pressure distribution data 
were introduced in a previous section. Section E-4 of this chapter mentionned the importance of a 
complete understanding of the panel numbering method in order to refer to the correct data. 


67 





68 















J W Q cn Q 
5 0 5 0 


30000000000000C 


3 O O O O O C 
3 O O O O O C 
3 O O O O O C 


w o 

300000000000000000000000 
300000000000000000000000"' 
300000000000000000000002XU 


ri000300000»-<rHrHr-li- 

?veooooi/>t^ 03 i-iou 3 r~o\r- 

niTiuisivoovoor-r-r^r-c^o 


"i<MCTvt^i<n'»a'n<-M^OOZO 

•^^oooooooooau 


< S 

a» O O 


5- m M D- O O c 

o «C u u II Z ti 

o z a H a: D C 
f-. J U} Q. o 3 
wwoe-ow- 
O O M b ec > 3 


o or-fiCTNOTfOr^r' rt^O'^ooooooiri rninooooor 

O 03 U3 O r-i <N r~ O n W) a <3^ ^ IN crn n 1*3 <M rH O 00 r 

O 00000000«-<>-l--<*-l»-<r'I'NrM<NfNrM/^rMfNrM<N(NfNIN»-*«- 

O - *00000000000000000000000000000 


ft. o 11 W ^ II X -11 • O o (h E X O O 11 

MU<OXOU iiUOCJii nf- 
mOi-iE-*XU-ujz II -< II oc >i Vi a. 0. lx. z t-f 
E-t-XM>XOi-<XXX<Q-o:t*>£i:o<J 

w«2e‘WMJxaxHmxouooo.m 

OJ3CZ(x:><ft.505oa:zzzzxz 


jO n O O O O O O O O O O O O O O O O O O O O O O O O O O o o 

-Z II MOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 
Ci-ixu-xoooooooooooooooooooooooooooo 

i. 5 e^j 2 . 

1 cn<i-ioooooooooooooooooooooooooooo 


o oo o o o 
O O o o o o 
o o O O O O I 


3000000000000000 


I_I_____ 

Figure 3.53 Bad PMARC Input File With Reverse Section Definition. F5 Wing, 
65A004.8 Section, AR=2.98. (Last Few Lines Not Shown). 
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G. 


TRIANGULAR WING 


The following wing is also called delta wing. We changed the terminology to triangular 
to ensure distinction from the delta wings studied earlier and also to match the literature which was 
used as reference. The triangular wing modeled here has a thickness distribution represented by the 
Equation (3.5). 

y . KxA ■ 26(x-M)(1.;c) ( 3 , 5 ) 

a 


which is a parabolic arc in all X,Y planes and which has linearly decreasing h in all Y,Z planes. 
(Note that the roll, pitch, and yaw axis correspond to the x, y, and z axis respectively in PMARC and 
to the X, z, and y axis respectively in our reference literature (Ref. [15]). Figure 3.54 illustrates the 
wing. Its aspect ratio is 2/3 and for the purposes of analytical and computer analysis, its root chord 
will be set at one. Using these data, it can be shown with basic relationships that the semi-span at 
the trailing edge is 1/6, and that it varies in the x-direction like x/6, where x is the non- 
dimensionalized root chord location. 



Figure 3.54 Triangular Wing, Parabolic Arc 
Section, AR=2/3. 
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Unlike the previous sections, the reference theory will be discussed in greater length before 
being compared to the PMARC results. The theory, introduced in the previous section, will first be 
presented in its general form, then we will analytically rework an example presented in the 
Reference in order to ascertain the validity and understanding of the theory. Once our results 
corroborate those of the author, we will finally proceed with the PMARC computation and 
comparison. As it will be seen, major disparities exist in the predicted velocity distributions. 

1. Theory, Bacl^round and Reference Data 

In the early fifties, Friedrich Keune (Ref. [15]) presented a theory on low aspect ratio wings 
with small thickness at zero lift in subsonic and supersonic flow. He showed that a first 
approximation of the flow around low aspect ratio wings can be obtained by dividing the flow into 
two components: the cross-sectional flow {phi) and the spatial influence {omega). 


= {(9(x,y^yQ{x,y^)]U^ (3.6) 


a. Approximate Theory 

The starting point of this approximate theory is the equation representing the exact 
velocity potential of a flow. With a thickness distribution over the wing plane represented by: 

y = ±h{x,z) (3.7) 

the source distribution in the wing area (Figure 3.55) is depicted by: 

. IhJX.OU’^ (3.8) 

dl 
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Therefore the exact velocity potential is described with the following double integral: 
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where F(x,z) is the wing area. 



Figure 3.55 Two-Dimensional Flow Around a Low-Aspect-Ratio Wing with 
Thickness but Zero Lift. (From Ref. [15]). 


When considering low aspect ratio wings, the ordinate zeta of the location of any 
source element may be considered small in comparison with the root chord of the wing, therefore 
allowing a Taylor series development of the exact equation for small zeta. Only the first order terms 
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are kept, the higher order terms being negligible. It is then possible to develop the expression into 
a two-dimensional (Equation (3.10)) and a rotational symmetric (Equation (3.11)) component of the 
potential flow . The derivation is detailed in Chapter H of Ref. [15]. 


— = - O(a^lna) 

[/«• nJsj(x) 


(3.10) 


^ r° h(x^) dz 

U°° An Jo 


(3.11) 


The two-dimensional potential flow (Equation (3.10)) is caused by a two- 
dimensional source distribution 2*Ui,rf*(d/dx)[h(x,zeta)] on a line of constant chord location and at 
a thickness location of zero, between two outer spanwise limits (Figure 3.56). These sources inflate 
the flat plate to the thickness y(z)=+-h(x=const,z), or to the cross-sectional area defining the section. 

The axisymmetric flow (Equation (3.11)) is induced by a source distribution along 
the x-axis within the wing (Figure 3.56). This flow is identical with that around a slender body of 
revolution. 


Now the sum of the two terms shown above represents the first order term of the 
approximation of the velocity potential in the incompressible flow: 
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Figure 3.56 Cross-Sectional Flow and Spatial Influence. (From Ref. [15]). 


The potential function for incompressible flow shown above is given for any 
distance from the wing. On the surface of the wing (y=0), the ordinate z has the same order of 
magnitude as the length zeta^ so that z must consequently be neglected with zeta in this 
approximation. This leads to a transformation of the potential function as detailed in Chapter IV of 
Ref. [15]. The velocity field of the cross-sectional flow at y=0 is then given by: 
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The u component in the above equation (phQ can also be expressed the following non-dimensional 
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The spatial influence is represened by: 


n 1 ^ ^ Nr/ / or M 1 r^o „ 
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The u component in the above equation (omega^) can also be written; 
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One marked particularity of this theory is that the spatial influence (axisymmetric 
flow) is only a function of the cross-sectional area. Therefore the wing can be represented by a 
simpler geometiy with equivalent sectional area (like its equivalent body of revolution) in order to 
simplify the calculations 

b. Analytical Application Example Using The Triangular Wing 

Using a computer analytical solver (Ref. [16]), the u component of both terms of 
the flow were separately computed for the triangular wing described at the beginning of this section. 
For this wing, the segment of the cross-sectional area for positive ordinates z is: 

P(x^) . 46o(1-x)(j:-M)J£L (3.24) 

2a a 


and the cross-sectional area is given by: 

Q(x) = 46ax\l-x) _ (3.25) 


(1) Spatial Influence. The u component of the spatial influence was 
computed using two equations. The first one is the general equation (Equation (3.23)) and the second 
one is the specific equation for the triangular wing (Equation (3.26)) as given by Ref. [15]. 


^ = -i~l- 4 -^(l-2x) - —(l-3x)ln4x(l-x) (3.26) 

6 a n l-x n n 
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Both equations would be expected to produce the same results as Equation (3.26) is directly derived 
from Equation (3.23). As Figure 3.57 shows however, small disparities are observed, especially past 
the mid-chord point. These differences are explained by the accuracy of the integration methods 
used. But more importantly, the figure also indicates that we were able to exactly reproduce the 
results of Ref. [15] using the general Equation (3.23). Appendix B presents the Equations (3.21) and 
(3.26) at different stages of their solution. 


Spatial Influence 



Figure 3.57 Comparison of Spatial Influence u-Component of The Velocity Potential as 
a Function of Chord Location. Triangular Wing, Parabolic Arc Spanwise Section, 
AR=2/3, AOA=0 Degree. 


(2) Cross-sectional flow. Once again, two equations were solved for 
comparison purposes. The specific equation (Equation (3.27)) as given in Ref. [15] and the general 
equation (Equation (3.17)). Appendix C presents the equations at different stages of their solution. 
The results were identical £ind in good agreement with those obtained by Ref. [15]. Figure 3.58 
presents these results. 
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Non-dimensionalized Spanwise Location (z/sigma) 


Figure 3.58 Comparison of The Chordwise Cross-Sectional Flow u-Component 
Velocity Potential as a Function of Non-Dimensionalized Spanwise Location. Triangular 
Wing, Parabolic Arc Section, AR=2/3, AOA=0 Degree. 
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(3) Total Flow. When the spatial influence and the cross-sectional flow 
results are added up, we obtain the curves depicted in Figure 3.59, which correspond very well to 
those computed by Ref. [15]. 


Total Flow 



Figure 3.59 Comparison of The Combined u-Components of The Chordwise Basic 
Flow as a Function of Non-Dimensionalized Spanwise Location. Triangular Wing, 
Parabolic Arc Spanwise Section, AR=2/3, AOA=0 Degree. 


2. Results and Discussion 

As for the F5 wing, PMARC seems not capable of recreating a realistic x-component 
velocity distribution on the triangular wing. Figure 3.60 once again shows velocity peaks that 
increase as the spanwise location gets towards the tip. This is a behavior that would be more 
representative of a wing of constant thickness. In our case, the tip is modelled to be of virtually zero 
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thickness, therefore the velocity perturbation should approach zero as we get closer to the tip. On 
Figure 3.60, zbar is the non-dimensional spanwise location divided by the dimensionless semi-span. 
As mentioned before, z corresponds to they axis in the PMARC notation and was left as is to match 
our reference notation (Ref. [15]). 



Figure 3.60 Spanwise Perturbation Velocity as a Function of Chord Location. 
Triangular Wing, Parabolic Arc Spanwise Section, AR=2/3, AOA=0 Degree. 


In order for the reader to better appreciate what behavior is expected, we extracted a figure 
from Ref. [15] which shows the perturbation velocity distribution on the triangular wing. 
Interpretation of this figure corroborates the intuitive expectation of progressively smaller 
perturbation velocity from root to tip. See Figure 3.61 below. 
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Figure 3.61 Isotaches in The Wing Plane. Triangular Wing, Parabolic Arc Spanwise 
Section, AR=2/3, AOA=0 Degree, z/x Proportions Not to Scale. (From Ref. [15]) 


To ensure that the erroneous results were not due to our geometry which in turn could have 
caused matrix singularities, a series of variation of our triangular wing were modeled. The tips were 
truncated, the leading and trailing edges were rounded, and the panel density and distribution were 
varied Larger aspect ratio delta wings were also modelled and all without any positive effect. The 
same erroneous behavior was always present. 

Interestingly however, the problem is not apparent in the case of rectangular geometries. 
As shown in Appendix D, a thickness tapered rectangular wing was modelled and the velocity 
perturbation decreases towards the tips, as it should. 
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3. 


PMARC Input Instructions 


The input file for the triangular wing is no different than any of the delta wings previously 
studied. Only the airfoil section, the tip location, the geometric fields, and the scaling vary. The 
points defining the airfoil cross-section were separately obtained substituting chordwise locations 
in the triangular wing expression of thickness distribution (Equation (3.5)). In the following 
paragraphs, one of the above differences will be discussed, the tip scaling, and an important 
observation will be brought to light regarding the spanwise cross-section definition in PMARC. 

a. The Tip Scaling 

As it can be observed in the triangular wing input file (Figure 3.62), the tip scaling 
was made smaller than what was done for the delta wings. It was mentioned then that we wanted as 
small a tip as we could get without creating too many panels with virtually zero area. In the case 
of the delta wings, SCALE=0.05 was an optimum value which would generate a relatively low 
number of zero area panels . With the triangular wing, this scaling would change the aspect ratio 
of the wing as it corresponds to a larger proportion of the root chord. We then had to go down as low 
as SCALE=0.01 to conserve the triangular shape, although that increased the number of zero area 
panels. However, these panels account for a small proportion of the total number of panels, still 
allowing for a good flow solution. 

b. Spanwise Cross-Section Definition 

The spanwise cross-section of the wing as defined in the input file of Figure 3.62 
has a linearly decreasing thickness. The cross-section would then have a diamond shape. This is due 
to the fact that we prescribed the coordinates for the origin of the root chordwise cross-section and 
also for the tip cross-section and from there, PMARC seems to linearly interpolate all the origin 
coordinates of all chordwise cross-sections in between. The delta wings studied earlier also had this 
diamond shape spanwise cross-section. 

If the spanwise cross-section has a non-linearly decreasing thickness from root to 
tip creating an oval shape or double parabolic arc shape, then the PMARC construction 


84 




























would be more intricate. Although not tried, it is suspected that several different airfoil section 
definitions would be required at various spanwise locations, each showing a different thickness to 
chord ratio. Enough of these intermediate section definitions would be required to create a smooth 
spanwise cross-section. 

4. PMARC Output - Data Retrieval 

a. Direct Method 

In order to reproduce the format of the results presented in our reference literature 
(Figure 3.63), an intricate data retrieval is required. As can be seen from that figure, we need arrays 
of velocities taken at constant jc values, for various increasing z values (y in PMARC). To do that, 
we need to find the beginning of the steady state results in DATA6. Then, for each block of results 
(which corresponds to columns, or strips of panels) the velocity Fx at x=0.1 is recorded. The same 
process is repeated for x=0.2 and all other chordwise locations desired. 

Unless the panel density is very high, one won't find in the output the chordwise 
locations that exactly correspond to locations plotted in Figure 3.63, (x=0.1, 0.2, 0.3, 0.4 ...). 
Therefore, an interpolation will be required. In that regard, it was observed that quadratic 
interpolations produce slightly different results than linear interpolations. The difference being 
small, the latter method would be acceptable for the type of comparative study in question. 

b. Indirect Method 

Another data retrieval technique is a lot simpler but some interpretation of the 
resulting plots is needed to compare with the reference plots. All that is required here is to plot the 
X column against the Vx column for selected blocks of data (or panel strips). This produces plots as 
presented in Figure 3.60. The difference is that instead of dealing with perturbation velocities as a 
function of spanwise stations for fixed chordwise locations, we deal with perturbation velocities as 
a function of chordwise stations for fixed spanwise locations. 
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Figure 3.63 Flow Results From F. Keune in The Plane y=0 (PMARC z). Triangular 
Wing, Parabolic Arc Spanwise Section, AR=2/3, AOA=0 Degree. (Fjom Ref. [15]). 


c. Perturbation Velocity 

The data retrieved in DATA6 are x-component velocities. Plotting these values 
instead of those shown in the reference plots (perturbation velocities) would only make a difference 
in the abscissa and ordinate; the behavior of the curves would be the same. However, we were 
interested in the perturbation velocities which are simply the velocities found in the DATA6 minus 
the freestream velocity, all divided by the product of the freestream velocity and the dimensionless 
maximum thickness. 


u 


Ub 

oo 


(3.4) 


87 






H. CONCLUSION 


In this chapter, we demonstrated that lift and pressure coefficients generated by PMARC 
accurately compare with the wide range of experimental and theoretical results consulted, or with 
higher order panel code results. However, disparities were found in the area of velocity distribution 
as the PMARC results showed increasing velocity perturbation towards the tips of the triangular 
wing whose thickness actually decreases towards the tips. On the other hand, it was also shown that 
the perturbation velocity behavior was adequate in the case of a thickness tapered rectangular wing. 
Another area that requires more attention is the wing tip modelling. It would have to be ascertain 
that panels with virtually no area indeed do not corrupt any of the output. 

The cases studied in this chapter include swept back wings, delta wings, rectangular wings 
and the F5 wing. Instructions were given on how to build each model and make the input file. Also, 
details on data retrieval from the PMARC output file were provided for each case studied. 
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IV. SOFTWARE VALIDATION - BODIES OF REVOLUTION 

A. INTRODUCTION 

The previous chapter studied the PMARC capability to recreate existing aerodynamic 
parameters of standard shape wings in steady state. This chapter has the same objective for bodies 
of revolution. The body of revolution is the longitudinal integration of circular cross-sections. 
Computational advantages are associated with this technique as was introduced earlier.The body of 
revolution can be one of independent shape or one which is an equivalent representation of the 
spanwise cross-sectional area of a given wing. The two types will be treated here. More specifically, 
the equivalent area body of revolution for the F5 wing and the triangular wing studied previously 
will be modelled. A spindle will also be analyzed. The spindle is a body of revolution symmetric 
about all axes originating at the mid-chord point. In all cases, the emphasis will be on the pressure 
distributions. It is because there exist important relationships between the aerodynamic coefficients 
of a wing and those of its equivalent body of revolution that we will build and analyze the latter. 
Note that the demonstration of the above mentioned relationships and their utilities will be the 
subject of follow-up work. For now, we will only present and discuss the. chordwise pressure 
distribution at various circumferential locations for angles of attack of zero and five degrees, which 
will be compared against the outcome of the slender body theoiy. 

Based on the format of the previous chapter, each body of revolution will have a dedicated 
section. The respective bodies will first be described, followed by an overview of the applicable 
theory serving as reference. Then, the results as obtained by PMARC will be presented and 
discussed. Also, the specific PMARC input steps for modeling the subject body of revolution will 
be explained. Finally, the output data retrieval techniques are presented. 
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B. F5 WING - EQUIVALENT AREA BODY OF REVOLUTION 


The F5 wing was modelled earlier in this paper and several pressure distributions were 
obtained for some specific conditions of flight. The correspondent equivalent area body of revolution 
will have the same main chord length of 0.6439 unit. Figure 4.1 shows the F5 wing superposed to 
its equivalent body of revolution. The spanwise cross-sectional areas of the wing needed to build the 
body of revolution were obtained using an indepth computer program designed by the author and 
presented later in this section. 



Figure 4.1 F5 Wing With Corresponding Equivalent Area 
Body of Revolution. 


1. Theory, Background and Reference Data 

In order to validate the PMARC results of pressure distribution for the F5 equivalent body 
of revolution, we will rely on past experience for the general expected behavior, and on the 
approximate slender body theoiy for the magnitude of the pressures. Our chordwise pressure 
distribution analysis will be performed at several circumferential locations to ensure consistency in 
the expected trends. 
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a. Expected General Behavior 

As a general behavior, the total lift on a body of revolution should be zero as the 
upward lift on the front portion of the body is cancelled by the downward lift on the rear portion. In 
terms of pressure distribution, it will be observed that for a change of angle of attack from zero to 
an arbitrary value, the pressure differential (from zero to a new value of angle of attack) at a specific 
aft chordwise location will be equal in magnitude but in the reverse direction than that of a related 
rearward location. This means that somewhere along the chord there is a point where this differential 
is bound to be zero. This point is not at the mid-chord, except for a symmetrical body of revolution 
like the spindle, as we will see later. 

The behavior described above would be modified for the chordwise distributions 
that lie along the 0 and 270 degree circumferential positions. Not only is the magnitude of the 
pressure differential for a zero to five degrees shift in angle of attack expected to be the same at two 
related forward and rearward positions, but also the direction of the pressure. 

b. Slender Body Theory 

According to the slender body theoiy, the pressure on a body at angle of attack is 
approximated by the pressure at zero angle of attack plus the pressure differential between its value 
at zero and at the new angle of attack. This differential can be calculated using the following slender 
body theory derivation in which R is the radius distribution and Q the cross-sectional area 
distribution. 


= ^Pa.O " ^Pda (41) 


= 2ai?'(x)cos6 (4.2) 
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where: 
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0(x) = 71.3x‘®-352.00x^743.17x*-875.60x’.630.85x^ 


-286.66X ^81.70x ^-14.13x^1.42x 2-0.04x 


(4.5) 


Using the expected behavior description given above and the numerical results of 
the slender body theory applied to the F5 wing equivalent body of revolution, we will now present 
the PMARC results. 

2. Results and Discussion 

Figures 4.3 to 4.8 presented below display the chordwise pressure distribution at various 
circumferential positions for two angles of attack, which will allow comparison with intuitive 
expectations . Figures 4.9 to 4.12 show the pressure perturbation when going from zero angle of 
attack to 5 degrees, as calculated by PMARC and the slender body theory. It should be noted that 
the circumferential positions of the chordwise body line of analysis were chosen to correspond with 
the control point coordinates of the chordwise strips of panels which were the closest to the points 
separating the body cross-section in 12 pie sections (Figure4.2). Finally, although the first and last 
point on all the figures were plotted, they should be disregarded. The reasons are not well 
understood but PMARC has difficulty with the pointy nose and tail when the model is built by 
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repetitive scaling like for the F5 and triangular wing equivalent bodies of revolution, by opposition 
to individual cross-section definition like for the spindle. These two methods will be reviewed later. 



Figure 4.2 Circumferential Angular Locations Studied on The 
F5 Wing Equivalent Body of Revolution Cross-Sections. 
(Dotted Lines) 


a. Chordwise Pressure Distributions 

The first figure (Figure 4.3) shows the expected results that the presswe distribution 
will be the same for any circumferential position when the body of revolution is at zero angle of 
attack. Figure 4.4 presents the same parameters but with the body at five degrees angle of attack. 
We now observe the progression in pressure from top to bottom at a particular aft section and 
regression at a particular rear section. This plot also shows the relative chordwise position of the 
minimum pressure for each chordwise distribution. 
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Figure 4.3 Chordwise Pressure Distribution for Various Circumferential Angles. F5 
Wing Equivalent Body of Revolution. AOA=0 Degree. 



Figure 4.4 Chordwise Pressure Distribution for Various Circumferential 
Angles. F5 Wing Equivalent Body of Revolution, AOA=5 Degrees. 
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Figures 4.5 to 4.8 compare the chordwise pressure distribution at various 
circumferential positions for zero and five degrees angle of attack. Each figure displays the 
variation for two corresponding circumferential positions (longitudinal mirror image) of the body. 
In all cases, the leeward pressure perturbation is equal and opposite to the windward pressure 
perturbation, as expected. At theta in the vicinity of 90 degrees, a pressure perturbation exists but 
remains approximatively constant along the chord. 




Figure 4.5 Comparison of Chordwise Pressure Distribution Along Selected 
Circumferential Angular Positions for Two Angles of Attack. F5 Wing Equivalent Body 
of Revolution. 
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Figure 4.6 Comparison of Chordwise Pressure Distribution Along Selected 
Circmnferential Angular Positions for Two Angles of Attack. F5 Wing 
Equivalent Body of Revolution. 




Figure 4.7 Comparison of Chordwise Pressure Distribution Along 
Selected Circumferential Angular Positions for Two Angles of Attack. F5 
Wing Equivalent Body of Revolution. 
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Figure 4.8 Comparison of Chordwise Pressure Distribution Along a Selected 
Circumferential Angular Position for Two Angles of Attack. F5 Wing Equivalent Body 
of Revolution. 


b. Pressure Perturbations 

Using the PMARC results on Figures 4.5 to 4.8 above, the difference in pressme 
between the body at zero and five degrees angle of attack was calculated for various circumferential 
positions and plotted against the pressure perturbation as computed with the slender body theory at 
the same respective circumferential positions and the same change of angle of attack (Figures 4.9 
to 4.12). 


For reasons given in the previous subsection, the extremities of these plots should 
be discarded. However, this restriction still leaves relatively large disparities between the two 
perturbation calculations for all the cases considered. These disparities are larger towards the nose 
and the tail and progressively reduce to zero towards the shifting point. 
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The shifting point is defined as the chordwise location where there is no pressure 
perturbation following the change of angle of attack. Ahead of the shifting point the pressure 
perturbations will be in the opposite direction than those after. It can be observed that the shifting 
points as calculated with the slender body theory do not always correspond to those obtained by 
PMARC. The differences for the theta angles 4.4,31.5, 58.5, 121.5, 148.5, and 175.6 degrees are 
1.5%, 0.0%, 12.7%, 3.1%, 3.3%, and 6.9% respectively. 

For the theta position of 86 degrees, we observe that the perturbations are close to 
zero. A reading at exactly 90 degrees would in fact reveal zero perturbation when computed with 
the slender body theory. 




X/C 

Figure 4.9 Comparison of Chordwise Pressure Perturbation at Selected 
Circumferential Angular Positions for a 0 to 5 Degree Change in AO A. F5 
Wing Equivalent Body of Revolution. 
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Figure 4.10 Comparison of Chordwise Pressure Perturbation at Selected 
Circumferential Angular Positions for a 0 to 5 Degree Change in AO A. FS.Wing 
Equivalent Body of Revolution. 
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Figure 4.11 Comparison of Chordwise Pressure Perturbation at Selected 
Circumferential Angular Positions for a 0 to 5 Degree Change in AO A. F5 Wing 
Equivalent Body of Revolution. 



Figure 4.12 Comparison of Chordwise Pressure Perturbation at Selected 
Circumferential Angular Positions for a 0 to 5 Degree Change in AOA. F5 Wing 
Equivalent Body of Revolution. 
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3. 


PMARC Input Instructions 


The basic steps to construct a body of revolution are very similar from those required for 
a wing. The main difference resides in the section definition of the circles distributed along the 
longitudinal axis. There exists more than one way to construct a body of revolution. One will be 
presented here and a second when introducing the spindle in a later section. A common requirement 
for each method is a set of spanwise cross-sectional areas sufficiently dense to properly define the 
body of revolution. Section 6 will present a computer program designed by the author that generates 
41 cross-sectional areas and the chordwise locations of their centroids from the geometry of the F5 
wing. Each of the two body-of-revolution construction methods presented have advantages and 
disadvantages which will be discussed at the end of each section. Figure 4.13 presents the input file 
for the first method and the next paragraphs discusses it. 

a. Basic Input Section 

The only fields which may differ from the standard wings studied so far are CBAR, 
SREF, and SSPAN. CBAR is taken to be the length of the longitudinal axis of the body or 
revolution. SREF is the reference area upon which all the coefficients are calculated, and therefore 
its selection could be different depending on the analysis performed or the results desired. In our 
case, it was the area of the largest cross-section. Finally, for the pressure and lift coefficient analysis 
we did up to now, the SSPAN field did not seem to be a useful factor. For completeness, we 
selected it to be half the diameter of the largest cross-section. 

b. Assembly and Component Coordinate System Information 

Nothing is different from the standard wings studied so far. 

c. Patch and Section Coordinate System Information 

This section is significantly different from what we have seen so far. This section 
also offers some flexibility which we will demonstrate by building the F5 wing body of revolution 
one way, and later the spindle using another way. 
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Oc{31 199S14:S8 


F5eqArea_DATA5 


Page 3 


fcSECTl 

STX* 0.6006, 

STY* 0.0000 


ALT* 0.0, 

THETA* 0.1 


INMODE* 0, 

TNODS* 0, 

tSECTl 

STX* 0.6028, 

STY* 0.0000 


ALT* 0.0, 

THETA* 0. 


INMODE* 0, 

TNODS* 0, 

&SECT1 

STX* 0.6057, 

STY* 0.0000 


ALT* 0.0, 

THETA* 0. 


INMODE* 0, 

TNODS* 0, 

SSECTl 

STX* 0.6086, 

STY* 0.0000 


ALF* 0.0, 

THETA* 0. 


INMODE* 0, 

TNODS- 0, 

fcSECTl 

STX* 0.6122, 

STY* 0.0000 


ALF- 0.0, 

THETA* 0. 


INMODE* 0, 

TNODS* 0, 

tSECTl 

STX* 0.6160, 

STY* 0.0000 


ALF- 0.0, 

THETA* 0. 


INMODE* 0, 

TNODS* 0, 

&SECT1 

STX* 0.6201, 

STY* 0.0000 


ALF* 0.0, 

THETA* 0. 


INMODE* 0, 

TNODS* 0, 

fcSECTl 

STX* 0.6247, 

STY* 0.0000 


ALF* 0.0, 

THETA* 0. 


ZNMODE* 0, 

TNODS* 0, 

tSECTl 

STX* 0.6291, 

STY* 0.0000 


ALT- 0.0, 

THETA* 0. 


INMODE* 0, 

TNODS* 0, 

tSECTl 

STX* 0.6341, 

STY* 0.0000 


ALF* 0.0, 

THETA* 0. 


INMODE- 0, 

TNODS* 0, 

SSECTl 

STX* 0.6389, 

STY* 0.0000 


ALF* 0.0, 

THETA* 0. 


INMODE* 0, 

TNODS- 0, 

tSECTl 

STX- 0.6439, 

STY* 0.0000 


ALF- 0.0, 

THETA* 0. 


INMODE- 0, 

TNODS* 5, 

&BFNODE 

TNODE-3, 

TNPC-0, T 

&WAKE1 

IDWAK*0, 

IFLXW-l, 

WING WAKE 


&WAXE2 

KWPACH-2, 

KWSIDE*4, 


KWPAN2-0, 

NODEW-0, 

&WAKE2 

KWPACH*!, 

KWSIDE*2 , 


KWPAN2-0, 

NODEW=5, 


SCALE= 0.0329, 


SCALE* 0.0299, 


SCALE* 0.0276, 


SCALE* 0.0252, 


SCALE* 0.0226, 


SCALE* 0.0195, 


SCALE* 0.0195, 


SCALE* 0.0138, 


SCALE* 0.0113, 


SCALE* 0.0080, 


SCALE* 0.0002, 


SCALE* 0.0000, 


KWLINE=0, 
INITIAL=0, 
KWLINE=0, 
INITIAL-0, 


fcONSTRM NONSL =0, KPSL = 


RN =37210000, 

vise 

= 0.00015'J 

^ NSLBL 

= 1/ 


&END 

NVOLR* 

0, 

NVOLC* 

0, 




&END 

XO* 

-0.1000, 

YO* 

1.5000, 

zo* 

-0.1000, 

INTVSR* 1, 

&SND 

XI* 

-0.1000, 

Yl= 

1,5000, 

21* 

0.1000, 

NPTl* 0, 

&END 

X2 = 

-0.1000, 

Y2 = 

1.5000, 

Z2* 

0.1000, 

NPT2* 20, 

&END 

X3* 

1.1000, 

Y3* 

1.5000, 

Z3* 

-0.1000, 

NPT3* 25, 

&END 

XRO* 

2.0000, 

YRO* 

2.0000, 

ZRO* 

0.0000, 

INTVSC* 1, 

&END 

XRl* 

4.0000, 

YRl* 

2.0000, 

2R1* 

0.0000, 



XR2* 

2.0000, 

YR2* 

2.0000, 

ZR2* 

1.0000, 


&END 

Rl= 

0.1000, 

R2 = 

1.0000, 

PHIl* 

0.0, 

PHI2*360.0, 

6 END 

NRAD- 

5, 

NPHI* 

12, 

NLEN* 

3, 


&END 

NSTLIN' 

*0, 






&END 

SXO* 

-2.0000, 

SYO* 

1.0000, 

szo* 

-0.5000, 



SU- 

0.0000, 

SD= 

10.0000, 

DS* 

0.2500, 

INTSL* 1, 

&END 

SXO- 

-2.0000, 

SYO* 

1.0000, 

SZO* 

-0.4000, 



SU* 

0.0000, 

SD* 

10.0000, 

DS* 

0.2500, 

INTSL* 1, 

&END 

SXO* 

-2.0000, 

SYO* 

1.0000, 

SZO* 

-0.3000, 



su* 

0.0000, 

SD* 

10.0000, 

DS* 

0.2500, 

INTSL* 1, 

&END 

SXO* 

-2.0000, 

SYO* 

1.0000, 

SZO* 

-0.2000, 



su* 

0.0000, 

SD* 

10.0000, 

DS* 

0.2500, 

INTSL* 1, 

&END 

SXO* 

-2.0000, 

SYO* 

1.0000, 

SZO* 

-0.1000, 



su* 

0.0000, 

SD* 

10.0000, 

DS* 

0.2500, 

INTSL- 1, 

SEND 


Figure 4.13 (Cont’d) PMARC Input File for The F5 Wing Equivalent Area Body of 
Revolution. 
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The patch fields are not varying from before. It is the section fields that do. In the 
first method, a circle of radius one is defined by 19 points. As shown on Figure 4.13, that circle is 
then scaled to correspond to the cross-section radius at a particular chordwise location of the body 
of revolution. Thus, we scale the circle to zero at the nose tip and progressively augment and then 
reduce the scaling back to zero along the longitudinal axis. This process is done by repeating the 
SECTl block as many times as we have cross-section definitions. 

The number of panels on the circumference is prescribed in the first break point 
node (BPNODE) following the initial circle definition. The value attributed to INPC with this 
method determines the number of panels on one half of the circumference, as the model is built by 
joining two patches (Fig 4.14). The number of panels longitudinally from nose to tail is prescribed 
in the last SECTl block by the TNPS value. The panel distribution is selected as before with the 
TINTC and TINTS fields next to the respective panel fields. 

One last note is that the TNODS value is set at 0 in all the SECTl blocks to indicate 
that they constitute the first or an intermediate section of the patch. SECTl of the last block is set 
to 5 to indicate the last section definition. 

d. Wake and Streamlines 

It is important for a body of revolution to turn the wake option off. That is done by 
setting IDWAK to 0 in the WAKEl block. Indeed, the body of revolution theoretically does not 
generate a wake. 
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Figure 4.14 PMARC Panels and Patches Definition for The Body of Revolution. 


4. Scaled Circles Model Construction Method Characteristics 

The main advantage of this method is its ease of programing. The input file is no longer than 
for a simple wing. A circle is defined once and it is scaled at will. PMARC will create two patches 
to model the body of revolution which are mirror images of each other. One drawback when 
compared to the next method is the apparent difficulty with the pointed nose and tail. Even with a 
low number of panels which allow a finite panel area at the nose and tail, the same erratic behavior 
is present in the output as the pressures at the nose and tail break down. The rest of the results 
however are good. 
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5. 


PMARC Output - Data Retrieval 


For the purpose of our study, only the pressure coefficient distributions at various 
circumferential positions were required. Since patch one and two generate the same results for the 
control points in mirror image positions, only the Cp from patch one were extracted. In order to 
obtain a Cp distribution all along the chord for the theta position selected, one reading per column 
of data is required. One or more columns of data can be skipped but care should be exercised as to 
conserve the proper smoothness of the distribution. In the DATA6 output file, one column of data 
details the coordinates and aerodynamic parameters of one circumferential strip of panels running 
from bottom to top for patch one, and from top to bottom for patch two, as shown on Figure 4.14. 

The circumferential position theta is translated to a panel number on the output file by 
looking through one column of data for the arctangent of the y over -z ratio which will give the 
angle needed. Once the position of the panel relative to the first panel of that column is determined, 
the corresponding panel in each column of data will be approximatively at the same angle theta. In 
our case, we wanted to study the pressure distribution at seven different circumferential positions: 
0, 30, 60,90,120, 150, and 180 degrees. As mentioned earlier, we did not have to study the mirror 
image position of these angles (0,330,300,270,240,210 and 180 degrees) as the results would be 
identical. Because of the automatic panel distribution option used in PMARC, there was virtually 
no way to make the control points correspond with the theta desired. Consequently, we had to select 
values close to the ones above: 4.4,31.5, 58.5, 86, 121.5, 148.5, and 175.6 degrees. But that did not 
matter for our analytical comparison as we used the same values in the slender body theory 
application. It should be noted however, that by varying the number of panels around the body, one 
would eventually find the divisions that best suit the requirements. 

6. Computing The Cross-Sectional Area of The Body of Revolution 

The following subsections will review the method employed to build the F5 wing equivalent 
body of revolution using the F5 wing PMARC output, which is our only source of data. We will first 
go through the steps indicating what data are required from DATA6, and how this selection can 
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become tremendously complicated for some types of wing, like the F5 wing, compared to simple 
rectangular wings. We will then discuss the program that was written to assist in the construction 
of the body of revolution, along with some of its output. 

a. Using The PMARC Output to Build The Body Of Revolution - Code 

The first objective while building an equivalent body of revolution is to obtain the 
section definition of the many spanwise cross-sections needed. Once the cross-sectional spanwise 
surface coordinates are identified for each 'cut', all that is left is to calculate their respective area. 
In the case of the rectangular wing, constructing spanwise cross-sections using the PMARC output 
is very easy and can be done in a matter of minutes. That is so because all the panel rows are 
parallel to the leading edge and the output provides the coordinates for each panel control point 
(Figure 4.15). It then becomes very simple to use the coordinates of the top and bottom row of 
panels to define the spanwise cross-section of this wing. In cases where the rows of panels are not 
parallel to the leading edge, like for the F5 wing, it becomes more complicated as the panel rows do 
not match with the spanwise cross-section anymore (Figure 4.15). 



Figure 4.15 PMARC Panel Distributions for The Rectangular and The F5 Wings. 


A series of steps is then required to create the data points for spanwise cross-sections fi-om the panel 
control points provided by the output file (DATA6) of PMARC. These steps are summarized below 
and each corresponds to a section of the MATLAB code written to digitally perform them. The 
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code is presented in Appendix E. It is to be noted that steps (1) to (5) were built to simplify the data 
retrieval required for this process. These first few steps reconstruct the chordwise airfoil cross- 
section along the span; an operation already performed by PMARC. The advantage of our method 
is that we only need to input the general non-dimensional airfoil cross-sectional definition once 
instead of having to input the several cross-sectional definitions required along the span in order 
to proceed with the determination of the spanwise cross-sections. 

(1) Chordwise Section Definition . The first step is to define the airfoil 
cross-section in three dimensions. The x and z coordinates (in the PMARC axis definition) are in 
this case taken from the input file (DATA5), themselves taken from the literature. They coordinates 
(spanwise locations of the various two-dimensional airfoil sections) is taken in the PMARC output 
file, DATA6. They correspond to the coordinates of the chordwise panel control points calculated 
by PMARC from the prescribed airfoil section and wing geomebty. Strictly on the basis of the above 
coordinates, all the following operations will be performed. 

(2) Chordwise Panel Width . Because the spanwise panels are not always 
the same width, we needed a routine that calculates the width of each chordwise strip of panels. This 
information is later used along with the sweep angle to determine the chordwise origin of the airfoil 
cross-sections covering the control points of each strip of panels. 

(3) Origin Of The Chordwise Panel Control Point Cross-Section . 
Trigonometric operations on the sweep angle and the width and cumulative width of the columns 
allow to determine the origin in global coordinates of each of the chordwise cross-sections that cover 
the control points of each strip of panels. 

(4) Non-Translated Coordinates Of The Control Points Cross-Sections . 
This step constructs fifteen coordinate sets of airfoil chordwise cross-section. Each cross-section 
is the representation of the corresponding panel strip control points. At this stage, all the sections 
constructed have their leading edge x value at zero. However, all y and z values are representative 
of the actual wing proportional geometry. 
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(5) Translating The Cross-Sections . Using the results of step c, each cross- 
section is now translated to its actual relative position, completing the chordwise cross-sectional 
definition of the wing. As mentioned earlier, these coordinates are available in DATA6, but the data 
management required to isolate the points needed is by far more tedious than using a single non- 
dimensional airfoil cross-section definition along with this program. 

(6) Fitting the Chordwise Sections Along The Span . In order to create 
spanwise cross-sections, we needed to have a means to identify the coordinates of any surface 
locations on a chordwise cross-section instead of only the coordinates of the panel control points. 
To achieve this goal, we fitted the top and bottom surface of each chordwise cross-sections with high 
order polynomials. This resulted in thirty-two polynomial expressions for one half wing. 

(7) Creating the Spanwise Cross-Sections . For different chordwise 
locations (forty-three in this case) and using the polynomials obtained above, we now determine 
different z coordinates for selected y locations (fifteen). This will create a series of data points 
defining forty-three spanwise cross-sections. Each of these spanwise cross-sections then was fitted 
with a polynomial in order to smooth out some small irregularities introduced by the previous curve 
fitting. 


(8) Plotting the Results. This step helps visualize the spanwise cross- 
sections from leading edge apex to trailing edge apex. Some of these plots are shown below. Figure 
4.16 shows eight of the forty-one spanwise cross-sections for one half wing. Figure 4.17 presents 
seven cross-sections that end at the wing tip parallel to the chord. The flat look of that tip reflects 
reality as this geometry allows for the mounting of a tip tank or missile. Figure 4.18 illustrates a 
selection of cross-sections located in the back portion of the wing past the tips. The curvature of 
some of these cross-sections is due to the very high order polynomial used in the program which 
were required to get convergence at the extremity between the top and bottom curves. These 
tortuous curves could in fact be replaced by straight lines. 
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Unit Vertical Location (z direction) 


0.02 


FWD Apex to Tip L.E 



Figure 4.16 Selected F5 Wing Spanwise Cross-Sections From Apex to Tip Leading 
Edge. (Half Wing). 
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0.01 s 


Tip L.E. to Tip T.E. 



Figure 4.17 Selected F5 Wing Spanwise Cross-Sections From Tip Leading Edge to Tip 
Trailing Edge. (Half Wing). 


Tip T.E. to Reer Apex 



Figure 4.18 Selected F5 Wing Spanwise Cross-Sections From Tip Trailing Edge to 
Rear Apex. (Half Wing). 
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(9) Calculating the Area of Each Spanwise Cross-Section. The final step 
before constructing the equivalent body of revolution is to obtain the area of each spanwise cross- 
section. Two methods were used in order to ascertain the results. The first one calculated the area 
between the curves representing the spanwise cross-sections (top and bottom surface respectively) 
and the abscissa by integrating the polynomials describing each curve. A translation of these curves 
was necessary before the integration was performed in order to eliminate any intersection with the 
abscissa which would have rendered erroneous areas. The second method used a trapezoid approach 
to directly get the area between the two curves defining each spanwise cross-sections. Both methods 
gave the same results, as expected. Appendix F summarizes these results. 

(10) Constructing the Fxjuivalent Body of Revolution . Now that all 
spanwise cross-sectional areas are known, the radius of the related equivalent area circles can be 
determined and located at the corresponding chordwise locations, thus producing the wing 
equivalent body of revolution (Figure 4.19). 



Figure 4.19 Wing Spanwise Cross-Section and Its Equivalent Area Circle. 
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(11) Cross-Sectional Area Distribution It is possible now to establish an 
expression describing the distribution of the spanwise cross-sectional area. All we need to do is to 
find the best fit for the arrays of data points defining the area as a function of chordwise locations. 
In the case of the F5 wing, this expression is the tenth order polynomial shown in Equation (4.5). 

b. Note on the Equivalent Body of Revolution Computer Code 

In the conceptual phase of the program described above, several data management 
complications arose that needed additional sub-routines that kept adding to its length and 
complexity. However, as this program was initially intended for a one-time use tool, no effort was 
expended in trying to keep it computationally efficient as the subroutines were getting added. The 
end result is a program that requires a relatively long running time, and experiences difficulty in 
plotting certain three-dimensional representations. 

This program is now thought to be an excellent tool for the application at hand. 
Providing some spanwise locations are prescribed as well as the sweep angle and the chordwise 
cross-sections scaling, this program will produce forty-three spanwise cross-sections. If it is to be 
used extensively in the future, it could be modified to operate more efficiently, therefore becoming 
quicker. The modification required is to combine all arrays of data produced by the code into a large 
matrix. All the variables would then hold double indices instead of one. This program could also 
work for any other types of wing for which the general non-dimensional chordwise cross-sectional 
coordinates are known. However, it would first have to be carefully reviewed in order to customize 
it to the wing under study. 

c. Camber of The Body of Revolution 

One aspect still needing some investigation is the question of the camber. If the 
wing camber needs to be reflected in the body of revolution, one possibility is to place the centroid 
of each spanwise cross-section equivalent area circles along the mean camber line of the wing. 
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c. 


TRIANGULAR WING - EQUIVALENT AREA BODY OF REVOLUTION 


The body of revolution making the subject of this section is made out of equivalent spanwise 
cross-sections of the triangular wing studied in Chapter Ed, section G. The longitudinal axis of the 
body of revolution will then correspond to the root chord of the wing which was made unity through 
non-dimensionalization. The spanwise cross-sectional area of the wing was a lot easier than for the 
F5 wing as an analytical expression for the cross-sectional area was given in Ref. [15] and 
reproduced in Equation (4.6) below. 


Q(x) = 46ax^(l-x) 


(4.6) 


Sigma is the dimensionless semi-span and is equal to 1/6. Delta is the dimensionless maximum 
thickness and will be set at 0.098. Figure 4.20 shows the PMARC representation of the body of 
revolution with its parent wing. 



Figure 4.20 PMARC Top and Side Representation of The Triangular Wing Equivalent 
Area Body of Revolution and its Parent Wing. Maximum Thickness. Ratio = 0.09. 

(Two Figures not at Same Scale). 
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1. 


Theory, Background and Reference Data 


In order to validate the PMARC results of pressure distribution for the triangular wing 
equivalent body of revolution, we will once again rely on past experience for the general expected 
behavior, and on the approximate slender body theory for the magnitude of the pressures. Like for 
the F5 equivalent body of revolution, our chordwise pressure distribution analysis will be performed 
at several circumferential locations to ensure consistency in the expected trends. 

a. Expected General Behavior 

As for the F5 equivalent body of revolution, the total lift on the triangular wing body 
of revolution should be zero as the upward lift on the fi-ont portion of the body is cancelled by the 
downward lift on the rear portion. In terms of pressure distribution, it will once again be observed 
that for a change of angle of attack from zero to 5 degrees, the pressure differential at a specific aft 
chordwise location will be equal in magnitude but in the reverse direction than that of a related 
rearward location. This means that somewhere along the chord there is a point where this differential 
is zero. This point is not at the mid-chord, except for a symmetrical body of revolution like the 
spindle, as we will see later. 

As for the F5 equivalent body of revolution, the behavior described above would 
be modified for the chordwise distributions that lie along the 90 and 270 degree circumferential 
positions. Not only is the magnitude of the pressure differential for a zero to five degrees shift in 
angle of attack expected to be the same at two related forward and rearward positions, but also the 
direction of the pressure. 

b. Slender Body Theory 

As already mentioned in the previous section, the slender body theory dictates that 
the pressure on a body at angle of attack is approximated by the pressure at zero angle of attack plus 
the pressure differential between its value at zero and at the new angle of attack. This differential 
can be calculated using the following slender body theory derivation in which R is the radius 
distribution and Q the cross-sectional area distribution. 
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(4.1) 


O’ - * Cp^, 


= 2a/?'(x)cos0 (4.2) 

R(x) = V(^) (4.3) 

71 

«'W - (4.4) 

271 71 

0(x) = 46ax2(l-x) (4.6) 


Using the expected behavior description given above and the numerical results of 
the slender body theory applied to the triangular wing equivalent body of revolution, we will now 
present the PMARC results. 

2 . Results and Discussion 

Figures 4.22 to 4.27 presented below display the chordwise pressure distribution at various 
circumferential positions for two angles of attack, which will allow comparison with intuitive 
expectations . Figures 4.28 to 4.31 show the pressure perturbation when going from zero angle of 
attack to a higher value, as calculated by PMARC and the slender body theory. It should be noted 
that the circumferential positions of the chordwise body line of analysis were chosen to correspond 
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with the control point coordinates of the chordwise strips of panels which were the closest to the 
points separating the body cross-section in 12 pie sections (Figure4.21). Finally, although the first 
point on all the figures was plotted, it should be disregarded. As mentioned earlier, the reasons are 
not well understood why PMARC has difficulty with the pointed nose and tail when the model is 
built by repetitive scaling like for the F5 and the triangular wing equivalent bodies of revolution, 
as opposed to individual cross-section definition like for the spindle. 



Figure 4.21 Circumferential Angular Locations Studied on The 
Triangular Wing Equivalent Body of Revolution Cross-Sections (Dotted 
Lines). 


a. Chordwise Pressure Distributions 

The first figure (Figure 4.22) shows the expected results that the pressure 
distribution will be the same for any circumferential position when the body of revolution is at zero 
angle of attack. Figure 4.23 presents the same parameters but with the body at five degrees angle 
of attack. We now observe the progression in pressure from top to bottom at a particular aft section 
and regression at a particular rear section. This plot also shows the relative chordwise position of 
the minimum pressure for each chordwise distribution. 
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Figure 4.22 Chordwise Pressure Distribution for Various Circumferential Angles. 
Triangular Wing Equivalent Body of Revolution, Delta=0.098, AOA=0 Degree. 
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Figure 4.23 Chordwise Pressure Distribution for Various Circumferential Angles. 
Triangular Wing Equivalent Body of Revolution, Delta=0.098, AOA=5 Degrees. 
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Figures 4.24 to 4.27 compare the chordwise pressure distribution at various 
circumferential positions for zero and five degrees angle of attack. Each figure displays the 
variation for two corresponding circumferential positions (longitudinal mirror image) of the body. 
In all cases, the leeward pressure perturbation is equal and opposite to the windward pressure 
perturbation, as expected. At theta in the vicinity of 90 degrees, a pressure perturbation exists but 
remains approximately constant along the chord. 



Figure 4.24 Comparison of Chordwise Pressure Distribution Along Selected 
Circumferential Angular Positions for Two Angles of Attack. Triangular Wing 
Equivalent Body of Revolution, Delta=0.098. 


119 













1 



1 


o.s 

o. 

o 


0 


0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

x/c 

Figure 4.25 Comparison of Chordwise Pressure Distribution Along Selected 
Circumferential Angular Positions for Two Angles of Attack. Triangular Wing 
Equivalent Body of Revolution, Delta=0.098. 
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Figure 4.26 Comparison of Chordwise Pressure Distribution Along 
Selected Circumferential Angular Positions for Two Angles of Attack. 
Triangular Wing Equivalent Body of Revolution, Delta=0.098. 
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Figure 4.27 Comparison of Chordwise Pressure Distribution Along a Selected 
Circumferential Angular Position for Two Angles of Attack. Triangular Wing Equivalent 
Body of Revolution, Delta=0.098. 


b. Pressure Perturbations 

Using the PMARC results on Figures 4.24 to 4.27 above, the difference in pressure 
between the body at zero and five degrees angle of attack was calculated for various circumferential 
positions and plotted against the pressure perturbation as computed with the slender body theory at 
the same respective circumferential positions and the same change of angle of attack (Figures 4.28 
to 4.31). 


For reasons given in the previous subsection, the nose results of these plots should 
be discarded. However, this restriction still leaves relatively large disparities between the two 
perturbation calculations for all the cases considered. These disparities are larger towards the nose 
and the tail and progressively reduce to zero towards the shifting point. 
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As described before, the shifting point is defined as the chordwise location where 
there is no pressure perturbation following the change of angle of attack. Ahead of the shifting point 
the pressure perturbations will be in the opposite direction than those after. It can be observed that 
the shifting points as calculated with the slender body theory do not always correspond to those 
obtained by PMARC. The differences for the theta angles 3.4, 32.6,62.6,117.4,147.4, and 176.6 
degrees are 8.5%, 0.0%, 225.0%, 17.0%, 0.0%, and 13.8% respectively. Unlike the F5 wing 
equivalent body of revolution where all the differences were within acceptable range, we now 
observe complete disagreement for the theta angle of 62.6 degrees. As before, the 30 degree range 
shows perfect agreement. For the theta position of 87 degrees, we observe that the perturbations are 
close to zero. A reading at exactly 90 degrees would in fact reveal zero perturbation when computed 
with the slender body theory. 




Figure 4.28 Comparison of Chordwise Pressure Perturbation at Selected 
Circumferential Angular Positions for a 0 to 5 Degree Change in AOA. Triangular Wing 
Equivalent Body of Revolution, Delta=0.098. 
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Figure 4.29 Comparison of Chordwise Pressure Perturbation at Selected 
Circumferential Angular Positions for a 0 to 5 Degree Change in AO A. 
Triangular Wing Equivalent Body of Revolution, Delta=0.098. 



Figure 4.30 Comparison of Chordwise Pressure Perturbation at Selected 
Circumferential Angular Positions for a 0 to 5 Degree Change in AO A. Triangular 
Wing Equivalent Body of Revolution, Delta=0.098. 
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Figure 4.31 Comparison of Chordwise Pressure Perturbation at a Selected 
Circumferential Angular Position for a 0 to 5 Degree Change in AOA. Triangular Wing 
Equivalent Body of Revolution, Delta=0.098 

3. PMARC Input Instructions 

The steps used to construct the PMARC model of the triangular wing equivalent body of 
revolution are the same as for the F5 wing equivalent body of revolution. The reader should refer 
to the previous section for the details. The only difference is in the acquisition of the spanwise cross- 
sectional areas. Unlike for the F5 wing, we did not need the assistance of a computer program to 
get them as an analytical expression was available (Equation (4.6)). Appendix G shows the results 
of the analytical expression at various chord location for use in the slender body theoiy. Figure 4.32 
presents the input file for the triangular wing equivalent body of revolution. 

4. Scaled Circles Model Construction Method Characteristics 

The same advantages and disadvantages of the method used for the construction of the F5 
wing equivalent body of revolution discussed in the previous section stand for the triangular wing 
equivalent body of revolution. 
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The data retrieval techniques used here to obtain the pressure coefficient distribution at 
different circumferential theta angles are the same as for the F5 wing body of revolution studied in 
the previous section. 
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D. THE SPINDLE 


The body of revolution making the subject of this section is symmetrical about the 
longitudinal axis and its mid-chord. The length of the longitudinal axis was made unity through non- 
dimensionalization. The thickness distribution represented by the radius distribution is shown in 
Equation (4.7) below. 


R(x) = 2x(x-x^) 


(4.7) 


where tau is the maximum thickness. Figure 4.33 shows the PMARC representation of the spindle. 



Figure 4.33 PMARC Representation of The Spindle. 


1. Theory, Background and Reference Data 

In order to validate the PMARC results of pressure distribution for the spindle, we will once 
again rely on past experience for the general expected behavior, and on the approximate slender 
body theory for the magnitude of the pressures. Like for the F5 and the triangular wing equivalent 
body of revolutions, our chordwise pressure distribution analysis will be performed at several 
circumferential locations to ensure consistency in the expected trends. 

a. Expected General Behavior 

As for the F5 and the triangular wing equivalent body of revolutions, the total lift 
on the spindle should be zero as the upward lift on the front portion of the body is cancelled by the 


127 







downward lift on the rear portion. In term of pressure distribution, it will once again be observed that 
for a change of angle of attack from zero to 5 degrees, the pressure differential at a specific aft 
chordwise location will be equal in magnitude but in the reverse direction than that of the 
corresponding rearward location (mirror image fi-om mid-chord). This means that this differential 
will be zero exactly at the mid-chord point, unlike for the F5 and triangular wing equivalent body 
of revolutions where it was further back. 

The behavior described above would be modified for the chordwise distributions 
that lie along the 90 and 270 degree circumferential positions. Not only is the magnitude of the 
pressure differential for a zero to five degrees shift in angle of attack expected to be the same at two 
related forward and rearward positions, but also the direction of the pressure. 

b. Slender Body Theory 

As already mentioned in the previous sections, the slender body theory dictates that 
the pressure on a body at angle of attack is approximated by the pressure at zero angle of attack plus 
the pressure differential between its value at zero and at the new angle of attack. Recall that this 
differential can be calculated using the following slender body theory derivation in which R is the 
radius distribution and tau the maximum thickness. 


Cp = Cp^.o ^ ^Pda (4-1) 

= 2aJ?I'(jc)cos0 (4.2) 

R(x) = 2x(x-x^) (4.8) 
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Jl'(x) = 2t(1-2x) 


(4.9) 


Using the expected behavior description given above and the numerical results of 
the slender body theory applied to the spindle, we will now present the PMARC results. 

2. Results and Discussion 

Figures 4.35 to 4.40 presented below display the chordwise pressure distribution at various 
circumferential positions for two angles of attack, which will allow comparison with intuitive 
expectations . Figures 4.41 to 4.44 show the pressure perturbation when going from zero angle of 
attack to a higher value, as calculated by PMARC and the slender body theory. 

It should be noted that the circumferential positions of the chordwise body line of analysis 
were chosen to exactly correspond with the control point coordinates of the chordwise strips of 
panels which separate the body cross-section in 12 pie sections (Figure4.34). 

This exact selection was not possible in the two previous cases. It was made possible here 
due to the method employed to build the model as will be discussed in the next subsection. Also as 
a result of that particular method, the end points (nose and tail) were not corrupted and generated 
valid results. 




Figure 4.34 Circumferential Angular Locations 
Studied on The Spindle Cross-Sections. 


a. Chordwise Pressure Distributions 

The first figure (Figure 4.35) shows the expected results that the pressure 
distribution will be the same for any circumferential position when the body of revolution is at zero 
angle of attack. Figure 4.36 presents the same parameters but with the body at five degrees angle 
of attack. We now observe the progression in pressure from top to bottom at a particular section aft 
of the shifting point and regression at a particular section rearward of the shifting point. T his plot 
also shows the relative chordwise position of the minimum pressure for each chordwise distribution. 

Figures 4.37 to 4.40 compare the chordwise pressure distribution at various 
circumferential positions for zero and five degrees angle of attack. Each figure displays the 
variation for two corresponding circumferential positions (longitudinal mirror image) of the body. 
In all cases, the leeward pressure perturbation is equal and opposite to the windward pressure 
perturbation, as expected. At 90 degrees, a pressure perturbation exist but remains approximately 
constant along the chord. 
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Figure 4.35 Chordwise Pressure Distribution for Various Circumferential Angles. 
Spindle, Delta=0.160, AOA=0 Degree. 
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Figure 4.36 Chordwise Pressure Distribution for Various Circumferential Angles. 
Spindle, Delta=0.160, AOA=5 Degrees. 




Figure 4.37 Comparison of Chordwise Pressure Distribution Along Selected 
Circumferential Angular Positions for Two Angles of Attack. Spindle, Delta=0.160. 


132 










Figure 4.38 Comparison of Chordwise Pressure Distribution Along Selected 
Circumferential Angular Positions for Two Angles of Attack. Spindle, Delta=0.160. 



Figure 4.39 Comparison of Chordwise Pressure Distribution Along Selected 
Circumferential Angular Positions for Two Angles of Attack. Spindle, Delta=0.160. 














Figure 4.40 Comparison of Chordwise Pressure Distribution Along a Selected 
Circumferential Angular Position for Two Angles of Attack. Spindle, Delta=0.160. 


b. Pressure Perturbations 

Using the PMARC results on Figures 4.37 to 4.40 above, the difference in pressure 
between the body at zero and five degrees angle of attack was calculated for various circumferential 
positions and plotted against the pressure perturbation as computed with the slender body theory at 
the same respective circumferential positions and the same change of angle of attack (Figures 4.41 
to 4.44). 


Like for the F5 and triangular wing equivalent body of revolution, there exist for 
the spindle some disparities between the two perturbation calculations for almost all the cases 
considered. These disparities are larger towards the nose and the tail and progressively reduce to 
zero towards the shifting point. 
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Once again, it can be observed that the shifting points as calculated with the slender 
body theory do not always correspond to those obtained by PMARC. The differences for the theta 
angles 0, 30, 60, 120, 150, and 180 degrees are 6.0%, 3.0%, 28.0%, 19.0%, 4.0%, and 8.0% 
respectively. Here, we observe the largest deviation at a theta of 60 degrees. It appears like the 
deviation gets larger as theta approaches 90 degrees from either side. 

For the theta position of 90 degrees, we observe that the perturbations are zero 
when computed with the slender body theory, but the PMARC results follow an elliptic distribution. 




x/c 


Figure 4.41 Comparison of Chordwise Pressure Perturbation at Selected 
Circumferential Angular Positions for a 0 to 5 Degree Change in AOA. Spindle, 
Maximum Thickness Ratio = 0.160. 
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Figure 4.42 Comparison of Chordwise Pressure Perturbation at Selected 
Circumferential Angular Positions for a 0 to 5 Degree Change in AOA. Spindle, 
Maximum Thickness Ratio = 0.160. 



Figure 4.43 Comparison of Chordwise Pressure Perturbation at Selected 
Circumferential Angular Positions for a 0 to 5 Degree Change in AOA. Spindle, 
Maximum Thickness Ratio = 0.160. 
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Figure 4.44 Comparison of Chordwise Pressure Perturbation at a Selected 
Circumferential Angixlar Position for a 0 to 5 Degree Change in AOA. Spindle, 
Maximum Thickness Ratio = 0.160. 


3. PMARC Input Instructions 

The steps used to construct the PMARC model of the spindle significantly vary from those 
used for the F5 and triangular wing equivalent body of revolution. We will present here the second 
method of modelling a body of revolution in PMARC. One similarity however is the acquisition of 
the cross-sectional areas. Although we did not need the assistance of the computer program as with 
the F5 wing, an analytical expression was available (Equation (4.7)). Figure 4.45 presents an 
extract of the input file for the spindle. The following subsections will review this second approach. 
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a. Basic Input Section 

No differences exist in the format and content of the basic input section of this 
method and the first one used for the F5 and triangular wing equivalent body of revolution. 

b. Assembly and Components System Information 

Once again, the two methods present no differences in this area of the input file. 

c. Patch and Section Coordinate System Information 

This section is the one that differs for the two methods. Instead of defining the 
coordinates of one circle and scaling the latter to the dimensions corresponding to the area at a 
prescribed chord location, this method redefines new circle coordinates at each of these chord 
locations. Like for the first method, the SECTl block is repeated each time a circle is defined. 
However, the reader will recall that in that method the field INMODE was set to zero in all 
succeeding blocks to indicate that the circle section defined in the first block was to be copied. It 
was then scaled by changing the SCALE field. 

In the second method, INMODE is set at four in all blocks to indicate that the x, y, 
z coordinates of each circle will be inputted. In this case the field SCALE is always set to zero since 
the circles are defined in global coordinates. 

In the approach presented above, all the panel fields (TNPS and TNPC) were set to 
zero. This tells PMARC to use the input section definitions to define the panel edges, longitudinally 
and along the circumference. This was possible here because all the circle points were generated 
so they longitudinally corresponded to each other. Otherwise, TNPS and TNPC can be prescribed 
and PMARC disregards the points coordinates. 

d. Wake and Streamlines 

As for the previous method, the wake shall be turned off for the spindle. As for the 
streamlines, they are not required. 
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4. Global Coordinate Circles Model Construction Method Characteristics 

For the studies performed in this paper, the only advantage that this method seems to present 
is that it allows PMARC to compute more easily the flow at the pointed ends. This is said with some 
reservation because until we understand why the behavior of the first method is encountered in the 
vicinity of the pointed ends, there is no certitude about the actual superiority of the second method. 
Maybe it is possible to correct the input file in the first method to allow for a smoother computation 
at the ends. That would make the first method better all around for the type of studies performed 
because the input file for the spindle method is about thirty times longer. Despite this apparent 
disadvantage, this approach could however be necessary in some other configurations. 

5. PMARC Output - Data Retrieval 

For the purpose of our study, only the pressure coefficient distributions at various 
circumferential positions were required. In order to obtain a Cp distribution all along the chord for 
the theta positions selected, one reading per column of data is required. One or more columns of 
data can be skipped but care should be exercised as to conserve the proper smoothness of the 
distribution. In the DATA6 output file, one column of data details the coordinates and aerodynamic 
parameters of one circumferential strip of panels running all around the body from bottom to bottom 
as shown on Figure 4.46. 



Figure 4.46 PMARC Panel Distribution on The Spindle Single Patch Model. 
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The circumferential position theta which need a pressure distribution analysis is translated 
to a panel number on the output file by looking through one column of data for the arctangent of the 
y over -z ratio which will give the angle needed. Once the position of the panel relative to the first 
panel of that column is determined, the corresponding panel in each column of data will be 
approximately at the same angle theta. In this case, we were able to study the pressure distribution 
at circumferential positions dividing the half body in six exact parts: 0, 30, 60, 90,120, 150, and 
180 degrees. As mentioned earlier, we did not have to study the mirror image position of these 
angles (0, 330, 300, 270,240, 210 and 180 degrees) as the results would be identical. Recall that 
it is the cross-section point definitions that were judiciously made so that we ended up with control 
points exactly at these theta angles. 

E. CONCLUSION 

From our study of three different bodies of revolution, it became apparent that PMARC 
generates pressure distributions that match the expected trends. When compared to the slender body 
theory, our results indicate good agreement near the chordwise location where the pressure 
differential between the zero and five degrees angle of attack switches sign. As we proceed towards 
either end of the body, the difference between PMARC results and the slender body theory in the 
computation of these pressure differentials (or perturbations) progressively gets larger to an 
unacceptable level at the ends. 

Interpretation of the plots presented also reveals that the shifting point is different depending 
on the method of calculation (PMARC or slender body theory). The differences however are in most 
cases in the vicinity of 5%. In all cases, PMARC is expected to be more precise than the slender 
body theory. 
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V. UNSTEADY CASE - OSCILLATION 


A. INTRODUCTION 

The intent for this chapter is to validate the PMARC results for wings and bodies of 
revolution in oscillation. Time constraint prevented a complete analysis of the unsteady case. It was 
possible to obtain results for the F5 and the triangular wings and their respective equivalent bodies 
of revolution as well as for the spindle for ten-degree amplitude oscillations about their mid-chord 
at various reduced frequencies. The results will be presented with a warning to the reader about their 
validity. It is left for future work to analytically or experimentally verify these results. 

B. F5 WING AND EQUIVALENT BODY OF REVOLUTION 

The F5 wing was shown and described in section F of Chapter HI, and its equivalent body 
of revolution in section B of Chapter IV. Recall that for the purpose of the PMARC modelling, the 
main chord was set at 0.6439 unit and the corresponding maximum thickness was 0.0309 unit. The 
dimensionless maximum thickness ratio is then 0.0480. The wing and its-equivalent body of 
revolution are reproduced below in Figure 5.1. 



Figure 5.1 F5 Wing and its Corresponding Equivalent Area Body of Revolution. 
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This section will first lead the reader through the various calculations of frequency, time¬ 
stepping and other parameters related to each oscillatory case. The results will then be presented 
with the lift and moment coefficients plotted for two cycles of oscillation. Following the results, the 
specific input instructions to run an oscillatory case in PMARC will be described. Finally, a 
technique will be shown to retrieve the numerous data required from DATA6. 

1. Oscillatory Frequency and Time-Stepping Calculations 

Three reduced frequencies (A:) were selected for the oscillations: 0.1, 0.5, and 1.0. The 
oscillations had an amplitude of plus and minus five degrees about the mid-chord. Two full 
oscillatoiy cycles were required. Here is a sample calculation for k=0.1 and the results for the other 
cases. All parameters shown below apply for both the wing and the equivalent body of revolution. 

In the following equations, c was taken to be the root chord, and the freestream velocity was 
arbitrarily selected at 11 units per second. Omega is the oscillatory frequency in rad/sec and its 
PMARC notation is WRY. T is the period, iVthe number of time-steps and dt the size of the time- 
step in seconds. The PMARC notations of the last two parameters are NTSTPS and DTSTEP 
respectively. 


yfc=0. \rad—= ^ - -^- ^(xi^WRYAjmradls 

U 11 

00 


(5.1) 


=3.67mec/cycle 

(D 1.7083 


(5.2) 


t=N*dt=NTSTFS*DTSTEP^3.67i0=NTSTPS*0.1^NTSTPS=36.7Ssteps/cycle (5.3) 


NTSTPS=73.56 time-steps for two cycles 
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As noticed, the number of time-steps (NTSTPS) times the size of each time-step (DTSTEP) 
must equal the period (T). The combination selected is the one which will render a smooth 
distribution. In PMARC, this becomes a trial and error exercise. The following table provides the 
values of the above parameters for all values of reduced frequency studied. 


k 

Amplitude 

Osc. About... 

Omega (WRY) 

NTSTPS 

DTSTEP 

0.1 

+- 5 Degrees 

Mid-Chord 

1.7083 rad/s 

74 

0.1 sec 

0.5 

+- 5 Degrees 

Mid-Chord 

8.5417 rad/s 

30 

0.05 sec 

1.0 

+- 5 Degrees 

Mid-Chord 

17.0834 rad/s 

37' 

0.02 sec 


Table 5.1 Test Parameters for The Oscillation of The F5 Wing Equivalent Body of Revolution. 


2. Results 

Figures 5.3 to 5.5 below show the combined CL and CM results of the wing and the 
equivalent body of revolution in oscillation at the three reduced frequencies listed in Table 5.1. The 
analysis of the relative maximums and minimums of the CL and CM curves as a function of the 
reduced frequency, as well as the general behavior will be left to future worL It will only be said 
that the lift and moment coefficients for the body of revolution seem to behave as expected. The 
lift vectors aft of the pivot point are balanced by those behind it (Figure 5.2 ). This causes the net 
value of the lift to be close to zero (zero in steady state) and the moments to increasingly become 
larger with angle of attack as larger forces each side of the pivot point and in opposite directions 
torque the body. Some lift is expected due to oscillation effects. All the oscillation figures are 
plotted against reduced time {(bar) which is obtained as follows: 


<0 = NTSTPS{i) * DTSTEP 






























Figure 5.2 Typical Balanced Pressure Distribution Along The 
Longitudinal Camber of a Body of Revolution. 




Figure 5.3 Comparison of The Lift and Moment Coefficients in Unsteady State as a 
Function of Reduced Time. F5 Wing and Equivalent Body of Revolution, Oscillation 
Amplitude=10 Degrees About c/2, Ujnf^ll, Omega=1.7083 rad/s, k=0.1. 
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Figure 5.4 Comparison of The Lift and Moment Coefficients in Unsteady State as a 
Function of Reduced Time. F5 Wing and Equivalent Body of Revolution, Oscillation 
Amplitude=10 Degrees About c/2, Uj„f=l 1, Omega=8.5417 rad/s, k=0.5. 
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Figure 5.5 Comparison of The Lift and Moment Coefficients in Unsteady State as a 
Fxmction of Reduced Time. F5 Wing and Equivalent Body of Revolution, Oscillation 
Amplitude=10 Degrees About c/2, Ui„^l 1, Omega=17.0834 rad/s, k=1.0. 
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3. 


PMARC Input Instructions 


In the unsteady case like in the steady case, the most difficult part about the conception of 
the input file is the geometry modelling. Throughout this paper we reviewed most of the techniques 
for simple wings and bodies of revolution. With this background, we got virtually all that it needed 
to run an unsteady case in PMARC. Only a few fields have to be changed to reflect the oscillation 
parameters. Theses fields are concentrated in the first two major sections of the input file: the Basic 
Input Section, and the Assembly and Component Coordinate Section. Figures 5.6 and 5.7 show 
these sections extracted from the input file of the F5 wing and its equivalent body of revolution 
respectively. The fields modified to run the oscillations are boxed. The following paragraphs will 
review these fields. The rest of the input file won't be discussed as there is no difference with the 
input file of the steady cases we ran previously. 

a. Basic Input Section 

In this section, four fields are directly related to the oscillatory case. They are: 
NTSTPS, DTSTEP, THEMAX, and WRY. Except for THEMAX, definitions and calculation 
methods for these fields were given in subsection 1 above. THEMAX is simply the amplitude of 
the oscillatory rotation about the PMARC global Y axis, in degree. Two other fields are not inherent 
to the oscillation but may need adjustment to allow a successful run. They are: MAXIT and 
SOLRES. The first one is the limit on number of solver iterations. The second is the convergence 
criteria for the matrix solver. The recommended setting are 150 to 200 for MAXIT and 0.0005 for 
SOLRES. 


b. Assembly and Component Coordinate System Information 
In this section of the input file, only one field requires to be changed from the steady 
case: ASEMX. That indicates in PMARC the axis of rotation for the oscillation. As can be noticed 
on Figures 5.6 and 5.7, this value needs to be negative in order to have the rotation axis on the body. 
That is because the body is defined in assembly coordinates and by putting a negative value to 
ASEMX we are prescribing the origin of the assembly coordinates to begin somewhere behind the 
global coordinates origin, the latter being taken for the axis of rotation. 


148 




























c. Customizing The Source Code If Failure to Run 
Two other parameters require attention. However, they are not found in the input 
file. They are in the parameter file of the source code (PARAM.DAT). One is the number of surface 
panels allowed (NSPDIM) and the other one is the number of wake panels allowed (NWPDIM). 
NSPDIM should at least be set to the number of panels covering the surface studied. NWPDIM was 
boosted up form 1900 to 5000 in order to get a successful run. Although it has not yet been 
confirmed or scientifically explained, it appears like this latter parameter should be set at a value at 
least equal to the product of the total number of spanwise panels and the number of time-steps. 
When the parameter file is modified, the source code has to be recompiled. 

4. PMARC Output - Data Retrieval 

In the oscillatory case, large files of output data are created; sometimes up to 6 Megabytes. 
Unless the number of time-steps is very small, the data retrieval should not be done manually as it 
would become very time consuming. The output file for the oscillation case follows the same format 
as for the steady case. At time-step zero and at the very last time-step, the x, y, z coordinates, the Cp, 
Vx, Vy, Vz, V, DUB, DPHI/DT and MACH number are provided for the control points of all panels. 
For the intermediate time-steps, the control point coordinates from the leading edge are given with 
the section chord, CL, CD, CY, Cm, Cn, and Cl. Figure 5.8 shows selected portions of a thirty page 
output file. 

Not shown on Figure 5.8 are the total coefficients for each time-steps. These are the values 
we needed in our study. More specifically CL and CM. In order to efficiently extract them, small 
programs written in the UNIX operating system AWK language were written. The data retrieval 
becomes relatively simple as the output file repeats the total coefficient block under the same format 
after each time-step. Thus, these programs look for a prescribed recurring feature, and from it 
counts an indicated number of lines and then take the desired value in the prescribed field of that 
line. In our case, the program was searching for the words WIND AXES from which it counted 32 
lines and then stored the first field for the CL value or the fourth for the CM value. This process was 
repeated at each time-step. These geometry specific programs can be found in Appendix H. 
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Figure 5.8 Selected Portion of a PMARC Unsteady Case Output File Showing The 
Data Produced at The Initial, First and Last Time-Step. 





















c. 


TRIANGULAR WING AND EQUIVALENT BODY OF REVOLUTION 


The triangular wing was shown and described in section G of Chapter m, and its equivalent 
body of revolution in section C of Chapter IV. For these cases, the main chord and longitudinal axis 
respectively were made unity. Like for the steady case, two thickness ratios will be treated in 
oscillation: 0.098, and 0.196. The wing and its equivalent body of revolution are reproduced below 
in Figure 5.9. Then the oscillation parameter calculations, the oscillation results, and the PMARC 
specifics are presented. 



Figure 5.9 PMARC Top and Side Representation of The Triangular Wing 
and Its Corresponding Equivalent Area Body of Revolution. Maximimi 
Thickness Ratio = 0.098. (Each View Not at Same Scale). 


1. Oscillatory Frequency and Time-Stepping Calculations 

The same three reduced frequencies (A:) as for the F5 were selected for these oscillation runs: 
0.1,0.5, and 1.0. The movement had the same plus and minus five degrees amplitude about the mid¬ 
chord, and two full oscillatory cycles were operated. Again, here is a sample parameters calculation 
for k=0.1 and the results for the other cases. All parameters shown below apply for both the wing 
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and the equivalent body of revolution, c is the root chord (or longitudinal axis), omega is the 
oscillatory frequency in rad/sec (WRY in PMARC), T is the period, N the number of time-steps and 
dt the size of the time-steps in seconds. The PMARC notations for the last two parameters are 
NTSTPS and DTSTEP respectively. The freestream velocity was arbitrarily selected at 11 unit per 
second. 


=. «=fFi?r=l.l rad/s 

U 11 

oa 


(5.4) 



=YY=5.71 SQcIcycle 


(5.5) 


T=N*dt=NTSTPS*DTSTEP - 5.7UNTSTPS*0.2 - NTSTPS=2i.55 stepsicycle 

(5.6) 


NTSTPS=57.1 time-steps for two cycles 


As mentioned in the previous section, the number of time-steps (NTSTPS) times the size 
of each time-step (DTSTEP) must equal the period (T). The combination selected is the one which 
will render a smooth distribution in the PMARC output. The following table now provides the value 
of the parameters shown above for all values of reduced frequency studied. Note that all the 
parameters not only apply for both the wing and the equivalent body of revolution, but also for the 
two thickness distributions considered. 


k 

Amplitude 

Osc. About... 

Omega (WRY) 

NTSTPS 

DTSTEP 

0.1 

+- 5 Degrees 

Mid-Chord 

1.1 rad/s 

58 

0.2 sec 

0.5 

+- 5 Degrees 

Mid-Chord 

5.5 rad/s 

23 

0.1 sec 

1.0 

+- 5 Degrees 

Mid-Chord 

11.0 rad/s 

23 

0.05 sec 


Table 5.2 Test Parameters for The Oscillation of The Triang. Wing Equivalent Body of Revolution. 
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2 . 


Results 


Figures 5.10 to 5.12 below show the combined CL and CM results of the wing and the 
equivalent body of revolution of thickness ratio 0.098 in oscillation at the three reduced frequencies 
listed in Table 5.2. Figures 5.13 to 5.15 do the same thing for a thickness ratio of 0.196. The 
analysis of the relative maximums and minimums of the CL and CM curves as a function of the 
reduced frequency, as well as the general behavior will be left to future work. It will only be said 
that the lift coefficients for the body of revolution seem to behave as expected only for the low 
reduced frequency, contrary to the F5 where their behavior seemed legitimate for all reduced 
frequencies. As illustrated in section B-2 of this chapter, the lift vectors in front of the pivot point 
are theoretically balanced by those behind it (Figure 5.2). This should cause the net value of the lift 
to be close to zero at any angle of attack. Some lift is expected due to the oscillation effects. Also, 
as observed in om results, the moments should increasingly become larger with angle of attack. As 
detailed in the previous section, all the figures are plotted against the reduced time which is the ratio 
of the product of freestream velocity and time over the main chord. 

3. PMARC Input Instructions 

Except for the values of the parameters, there is no difference between the input files for the 
oscillations of the triangular wing and its equivalent body of revolution and those of the F5 wing 
and its equivalent body of revolution studied in the previous section. One input file for the triangular 
wing and one for its equivalent body of revolution for the reduced frequency of 0.1 are presented 
below. The reader should consult section B-3 of this chapter for details. 

4. PMARC Output- Data Retrieval 

As no difference exists with the data retrieval techniques used in the previous section, the 
reader should refer to the latter for details on the program required to extract the data off DATA6. 
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Figure 5.10 Comparison of The Lift and Moment Coefficients in Unsteady State as a 
Frmction of Reduced Time. Triangular Wing of AR=2/3 and Equivalent Body of 
Revolution of Delta=0.098, Oscillation Amplitude=10 Degrees About c/2, Ui„f=l 1, 
Omega=l.l rad/s, k=0.1. 





Reduced Time 

Figure 5.11 Comparison of The Lift and Moment Coefficients in Unsteady State as a 
Fimction of Reduced Time. Triangular Wing of AR=2/3 and Equivalent Body of 
Revolution of Delta=0.098, Oscillation Amplitude=10 Degrees About c/2, Ui„f=l 1, 
Omega=5.5 rad/s, k=0.5. 
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Figure 5.12 Comparison of The Lift and Moment Coefficients in Unsteady State as a 
Function of Reduced Time. Triangular Wing of AR=2/3 and Equivalent Body of 
Revolution of Delta=0.098, Oscillation Amplitude=10 Degrees About c/2, Ui„f=l 1, 
Omega=l 1.0 rad/s, k=1.0. 
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Figure 5.13 Comparison of The Lift and Moment Coefficients in Unsteady State as a 
Function of Reduced Time. Triangular Wing of AR=2/3 and Equivalent Body of 
Revolution of Delta=0.196, Oscillation Amplitude=10 Degrees About c/2, Ujnf^l 1, 
Omega=l.l rad/s, k=0.1. 
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Figure 5.14 Comparison of The Lift and Moment Coefficients in Unsteady State as a 
Fimction of Reduced Time. Triangular Wing of AR=2/3 and Equivalent Body of 
Revolution of Delta=0.196, Oscillation Amplitude=l 0 Degrees About c/2, Uinf=l 1, 
Omega=5.5 rad/s, k=0.5. 
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Figure 5.15 Comparison of The Lift and Moment Coefficients in Unsteady State as a 
Fimction of Reduced Time. Triangular Wing of AR=2/3 and Equivalent Body of 
Revolution of Delta=0.196, Oscillation Amplitude=10 Degrees About c/2, Ui„f=l 1, 
Omega=l 1.0 rad/s, k=1.0. 
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Figure 5.17 PMARC Input File for The Triangular Wing Equivalent Body of 
Revolution in Oscillation About The Mid-Chord for k=0.1. 
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tSECTl 

STX* 0.7750, 

STY- 0.0000, 

STZ- 0.0000, 

SCALE* 0.0530, 



ALF* 0.0, 

THETA* 0.0 

t 




XNMODE- 0, 

TNODS= 0, 

TNPS* 0, TINTS* 

0, 

&END 

tSECTl 

STX= 0.8000, 

STY* 0.0000, 

STZ* 0.0000, 

SCALE* 0.0516, 



ALF= 0.0, 

THETA* 0.0 





INMODE= 0, 

TNODS* 0, 

TNPS* 0, TINTS* 

0, 

fcEND 

fiSECTl 

STX= 0.8250, 

STY* 0.0000, 

STZ* 0.0000, 

SCALE* 0.0498, 



ALF= 0.0, 

THETA* 0.0 





INMODE= 0, 

TNODS* 0, 

TNPS* 0, TINTS* 

0 t 

fcEND 

iSECTl 

STX= 0.8500, 

STY* 0.0000, 

STZ* 0.0000, 

SCALE* 0.0475, 



ALF= 0.0, 

THETA* 0.0 





INMODE= 0, 

TNODS* 0, 

TNPS* 0, TINTS* 

0/ 

fcEND 

&SECT1 

STX= 0.8750, 

STY* 0.0000, 

STZ* 0.0000, 

SCALE* 0.0446, 



ALF= 0.0, 

THETA* 0.0 





INMODE== 0, 

TNODS* 0, 

TNPS* 0, TINTS* 

0, 

fcEND 

&SECT1 

STXs= 0.9000, 

STY* 0.0000, 

STZ* 0.0000, 

SCALE* 0.0410, 



ALF= 0.0, 

THETA* 0.0, 





INMODE= 0, 

TNODS* 0, 

TNPS* 0, TINTS* 

0, 

&END 

&SECT1 

STX* 0.9250, 

STY* 0.0000, 

STZ* 0.0000, 

SCALE* 0.0365, 



AU= 0.0, 

THETA* 0.0, 





INMODE= 0, 

TNODS* 0, 

TNPS* 0, TINTS* 

0/ 

fcEND 

&SECT1 

STX= 0.9500, 

STY* 0.0000, 

STZ* 0.0000, 

SCALE* 0.0306, 



ALF- 0.0, 

THETA* 0.0, 





INMODE= 0, 

TNODS* 0, 

TNPS* 0, TINTS* 

0/ 

fcEND 

&SECT1 

STX= 0,9625, 

STY- 0.0000, 

STZ* 0.0000, 

SCALE* 0.0269, 



AIJ?‘= 0.0, 

THETA* 0.0, 





INMODE= 0, 

TNODS- 0, 

TNPS* 0, TINTS* 

0, 

fcEND 

tSECTI 

STX= 0.9750, 

STY- 0.0000, 

STZ* 0.0000, 

SCALE* 0.0222, 



ALr= 0.0, 

THETA* 0.0, 





INMODE= 0, 

TNODS- 0, 

TNPS* 0, TINTS- 

0, 

fcEND 

tSECTl 

STXs= 0.9813, 

STY- 0.0000, 

STZ* 0.0000, 

SCALE- 0.0194, 



ALF* 0.0, 

THETA- 0.0, 





INMODE= 0, 

TNODS* 0, 

TNPS* 0, TINTS* 

0, 

SEND 

fiSECTl 

STX* 0.9876, 

STY- 0.0000, 

STZ- 0.0000, 

SCALE- 0.0159, 



ALT* 0.0, 

THETA* 0.0, 





INMODE* 0, 

TNODS* 0, 

TNPS* 0, TINTS* 

0, 

fcEND 

tSECTl 

STX* 0.9938, 

STY- 0.0000, 

STZ* 0.0000, 

SCALE* 0.0113, 



ALF* 0.0, 

THETA* 0.0, 





INMODE* 0, 

TNODS* 0, 

TNPS* 0, TINTS* 

0, 

fcEND 

&SECT1 

STX* 0.9969, 

STY* 0.0000, 

STZ* 0.0000, 

SCALE* 0.0080, 



ALF* 0.0, 

THETA* 0.0, 





INMODE* 0, 

TNODS* 0, 

TNPS* 0, TINTS* 

Or 

6 END 

fcSECTl 

STX* 1.0000, 

STY* 0.0000, 

STZ* 0.0000, 

SCALE* 0.0000, 



ALF* 0,0, 

THETA* 0.0, 





INMODE* 0, 

TNODS* 5, 

TNPS* 22, TINTS= 

= 2, 

fcEND 

&BPNODE 

TNODE*3, 

TNPC=0, TINTC*0, SEND 




&WAKE1 IDWAK-0, irLXW=l, ITRFTZ- 0, INTRW= 0, ^ fcEND 

WING WAKE 


4WAKE2 KWPACH=2, KWSIDE^4, KWLINE=0, KWPAN1=0, 

KWPAN2=0, NODEW=0, INITIAL=0, SEND 


&WAKE2 

KWPACH*!, 

KWSIDE*2, 

KWLINE*0, 

KWPANl-0, 



KWPAN2*0, 

NODEW-5, 

INITIAL*0, 


&END 


iONSTRM 

NONSL 

=0, KPSL 

- 





fcEND 

fcBLPARAM 

RN *35755000, 

Vise 

= 0.cools'? 

, NSLBL = 1, 


&END 

&VS1 

NVOLR* 

0, 

NVOLC* 

0, 




fcEND 

&VS2 

X0= 

-0.1000, 

Y0 = 

1.5000, 

ZO* 

-0.1000, 

INTVSR* 1, 

&END 

&VS3 

XI* 

-0.1000, 

Yl= 

1.5000, 

Zl* 

0.1000, 

NPTl* 0, 

&END 

fcVS4 

X2 = 

-0.1000, 

Y2* 

1.5000, 

Z2= . 

0.1000, 

NPT2* 20, 

(END 

&VS5 

X3* 

1.1000, 

Y3* 

1.5000, 

Z3* 

-0.1000, 

NPT3* 25, 

&END 

fcVS6 

XRO* 

2.0000, 

YRO* 

2.0000, 

ZRO* 

0.0000, 

INTVSC* 1, 

&END 

&VS7 

XRl* 

4.0000, 

YRl* 

2.0000, 

ZRl* 

0.0000, 




XR2* 

2.0000, 

YR2* 

2.0000, 

ZR2* 

1.0000, 


fcEND 

&VS8 

Rl* 

0.1000, 

R2 = 

1.0000, 

PHIl* 

0.0, 

PHI2-360.0, 

fcEND 

fcVS9 

NRAD* 

5, 

NPHI* 

12, 

NLEN* 

3, 


&END 

fcSLINl 

NSTLIN^ 

=0, 






SEND 

fcSLIN2 

SXO* 

-2.0000, 

SYO* 

o 

o 

o 

o 

SZO* 

-0.5000, 




Figure 5.17 Cont'd PMARC Input File for The Triangular Wing 
Equivalent Body of Revolution in Oscillation About The Mid-Chord for 
k=0.1. 
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D. SPINDLE 


The spindle was shown and described in section D of Chapter IV. The main chord (or 
longitudinal axis) was made unity. The thickness ratio considered is 0.160. The spindle is 
reproduced below in Figure 5.18. Then the oscillation parameter calculations, the oscillation 
results, and the PMARC specifics are presented. 



Figure 5.18 PMARC Representation of The Spindle. 


1. Oscillatory Frequency and Time-Stepping Calculations 

This time, only two reduced frequencies (A:) will be studied: 0.1 and 0.5. The oscillation had 
the same plus and minus five degrees amplitude about the mid-chord, and two full oscillatory cycles 
were also computed. Again, here is a sample parameters calculation for k=0.1 and the results for 
the other case, c is the root chord (or longitudinal axis), omega is the oscillatory frequency in 
rad/sec (WRY in PMARC), T is the period, N the number of time-steps and dt the size of the time- 
steps in seconds. The PMARC notations for the last two parameters are NTSTPS and DTSTEP 
respectively. The ffeestream velocity was arbitrarily selected at 100 unit per second for the case of 
k=0.1 and at 11 for the case of k=0.5. 
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=» a)=priZy=10.0 rcuHs 


(5.7) 


U 100 


T~~=0.6283 secicycle (5.8) 

T=N*dt=NTSTPS*DTSTEP =» 0.6283=ATOrP5*0.025 =* NTSTPS=25A3 steps/cycle 

(5.9) 


NTSTPS=50.27 time-steps for two cycles 

As mentioned before, the number of time-steps (NTSTPS) times the size of each time-step 
(DTSTEP) must equal the period (T). The combination selected is the one which will render a 
smooth distribution in the PMARC output. The following table now provides the value of the 
parameters shown above for the two values of reduced frequency studied. Recall that the freestream 
velocity was 100 unit per second for A:=0.1 and 11 unit per second for ^0.5. This selection was 
necessary in order to keep the time step size from getting too small. 


k 

Amplitude 

Osc. About... 

Omega (WRY) 

NTSTPS 

DTSTEP 

0.1 

+- 5 Degrees 

Mid-Chord 

10.0 rad/s 

51 

0.025 sec 

0.5 

+- 5 Degrees 

Mid-Chord 

5.5 rad/s 

115 

0.02 sec 


Table 5.3 Test Parameters for The Oscillation of The Spindle. 


2. Results 

Figures 5.19 and 5.20 below show the combined CL and CM results for the two reduced 
frequencies considered. The analysis of the relative maximums and minimums of the CL and CM 
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curves as a function of the reduced frequency, as well as the general behavior will be left to future 
work. It can be said however that the behavior of the lift and moment coefficients warrants 
interrogation as to their validity. As illustrated in the section B-2 of this chapter, the lift vectors in 
front of the pivot point are theoretically balanced by those behind it (Figure 5.2). This should cause 
the net value of the lift to be close to zero at any angle of attack and the moment to vary over wider 
ranges. Some lift is expected due to the oscillation effects. Although the scaling of the plot varies 
widely, we still observe some lift, but more importantly, an equal moment variation. As detailed 
in section B of this chapter, all the oscillation figures are plotted against the reduced time which is 
the ratio of the product of the freestream velocity and time over the main chord. 



Figure 5.19 Lift and Moment Coefficient in Unsteady State as a Function of Reduced 
Time. Spindle, Delta=0.160, Oscillation Amplitude = 10 Degrees About c/2, Uj^lOO, 
Omega=10.0 rad/s, k=0.1. 
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Figure 5.20 Lift and Moment Coefficient in Unsteady State as a Function of Reduced 
Time. Spindle, Delta=0.160, Oscillation Amplitude = 10 Degrees About c/2, Ui„f=l 1, 
Omega=5.5 rad/s, k=0.5. 


3. PMARC Input Instructions 

Except for the values of the parameters, there is no difference between the input files for the 
oscillations of the spindle and those of the previous cases. An extract from one input file for the 
spindle is presented below. The reader should consult section B-3 of this chapter for further details. 

4. PMARC Output - Data Retrieval 

As no difference exists with the data retrieval techniques used in the previous oscillation 
sections, the reader should refer to the latter for details on the program required to extract the data 
offDATA6. 


165 




o o a a o c: 
z X z X X a 
u u M u «a » 


Q O Q 

z z z 

u u u 

«Jt •« 


SI4 


m 

>4 (II «a 

h) 


M I 

(0 Q 

< O 


o 00 
o 00 
• %o o 

« -O (f> • • 

00 00 

001* 

• • M u R > 

O O & M M N 
I I X Q M H 
N M O Q 0 < X 
N M U X U U 


^ O 


%% 


o o o o o o c 


300000000000000000000 




HXoooin -uo 

2 2 I R . 00 o Z « 
i N X O r- I O 
hi-iz 2 r •nmu 
x(nzxNobuz 
p.a. gg |zg 


o ~ 
0000 
o .00c 
« 0000 

a . * . o - 
a o o o • . 


a I o o o o 
. Q ■ ^ . • 
a Z O R o O c 



I o 

>4 I 
X H 

u u 

to X 

< H 


m 

S5“ 


R R k 
' H R R k 

X X ^ a 

^ H 04 X 

u u u o 


c « « « tn m ^ 


N O 00 "O ^ 


o*»'^irtin««'o'ov 


oooooooooooooooooooooooooooooooooooooooo 


^ o 000 


iT 

_ 


O O - "0 

o o o o C 

o o o o ► 

« -o « o o 
00 .^00 
000 • • c 

. . 00 

O O R R 

1 t X i 4 I a I 

X X o, 2 X X a 

X X r o X X S 

a, X o (o a< X D 

2 3 u o o o h 


x: o 
o o 
*i o 

« o 

ao 


< (0 r 
M Q «• 
b 3 O o 


0000 

UUUkRMUUkI 

ovw^to^inokir 

cM-x<eoe 9 win«> 

o 6 r*(o»r»r»^- 


00000000 


>00 
1 I I 
] ta u 


>000000000000000000 


I I I I I I I I I I 


' 4 <MO«' 0 *no\"i 




>000000000 


O (O *N (f) 

«» \0 r> «o o Ok 01 

0 Ok ko kfi 


klRlHUUMXUkl 

^ 0 k 040 k 0 «k 00 ku> 

k/>r«i*-i' 00 ®«U>*n 

r'Ok®<NOOr>itn® 

( 4 io®( 0 '*> 0 'weo-( 


U U kl U fc 

in O O m (/ 
«nwo«nr 
MraMfi*- 

o ^4 >4 o a 
kO <0 ko w << 


00000000000000000000000000000 


- 2 


S S ^ 


3 69 M 
COO 
coo 
H O O 

o o 
o o 
o o 


U M M «} 
0000 
0000 
0000 
0000 
0000 
0000 


>000 
>000 
^ ^ 


000 

000 

000 


0000 c 

0000 c 

-f 4- -*• + -I 

u u 
o o 
o o 
o o 
o o 


kr A n o n 03 0 


000000000 


> o o o o o c 

>000000000000000000000 

>000000000000000000000 

-—0000000000000000000 

ooooooooooooooooooo 
0000000000000000000 
ooooooooooooooooooo 


00c 
00c 
o o c 
0000 


bJ A A » 
0000 
0000 
0000 
0000 
0000 
0000 


00000000000000000000000000000 


Figure 5.21 Extract From The PMARC Input File for a Spindle in Oscillation About 
The Mid-Chord, Amplitude=10 Degrees, 0mega=10.0 rad/s, k=0.1. 
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E. CONCLUSION 


In this section, the lift and moment coefficients variations at different reduced frequencies 
were presented for the F5 wing, the triangular wing, their respective equivalent body of revolution, 
and the spindle. The magnitude and relative behavior of the results seem legitimate in some cases 
but require more study in others. 

In the event that the results are indeed erratic, one should not exclude the software as a 
potential source of error. However, the input files could also be introducing unrealistic behavior. 
One should then look for fields not yet discussed, or in the source code and its parameter sub¬ 
routines. 
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VI. SUMMARY AND CONCLUSION 


It was the objective of this thesis to contribute to the development of a new approach to 
determine the steady and unsteady airloads on missiles and on low-aspect ratio aircraft 
configurations typically used for high-performance fighter aircraft. This approach uses a second- 
order approximation to the steady and oscillatory lifting flow past low aspect ratio wings which 
exhibits the same fundamental division as the slender-body/slender-wing theory: incompressible 
cross-flow and a Mach-number dependent spatial influence. To test the range of validity of this 
second-order theoiy, it was necessary to compare it with an exact three-dimensional flow solution, 
and a low speed three-dimensional flow panel code named PMARC developed by NASA-Ames was 
selected for that goal. 

Therefore, PMARC was first validated by comparing it with other known solutions. For 
this purpose lifting flow over swept back, rectangular, and delta wings was studied and the PMARC 
predictions were found to be in good agreement. A detailed analysis of the F5 wing was then 
performed, and a program was written to compute its equivalent area body of revolution. The 
lifting pressure distributions were again found to agree well with the measufed ones. However, a 
discrepancy was detected where the PMARC predictions of the non-lifting pressure distributions on 
delta wings at zero angle of attack were compared with linearized theory and with the Oswatitsch- 
Keune low aspect ratio wing theoiy. PMARC predicts velocity perturbations which do not decrease 
with decreasing local chordwise thickness ratio of the wing. Interestingly, no similar discrepancy 
was found on lifting and non-lifting rectangular wings when thickness ratio decreases toward the 
tips. Similarly, the PMARC predictions for lifting and non-lifting bodies of revolution generally 
compare well with slender body theory. Some differences in the magnitude of the pressure 
perturbation exist towards the tips of the bodies of revolution, which are in some cases larger than 
10 percent. The chordwise location where the pressure perturbation switches sign also varies from 
the PMARC to the slender body theory, but never by more than 5 percent. Finally, PMARC was 
applied to the computation of the forces and moments generated on oscillatory wings and bodies of 
revolution. However, time constraint prevented any analytical or experimental comparison. 
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In light of the above summary, our contribution to the development of a second-order 
approximation to the steady and oscillatory lifting flow past low aspect ratio wings may appear dim 
as we were not able to obtain correct velocity distributions around delta wings at zero angle of 
attack using PMARC. However, through this difficulty we fulfilled one of our objectives which 
was to validate the software. It proved good in all other areas. It remains to work towards correcting 
PMARC or the procedure that lead to the erroneous results. Otherwise it will be necessary to 
establish the validity of another software. 


170 



APPENDIX A. F5 WING SPATIAL INFLUENCE FLOW CALCULATION 


This appendix shows all the intermediate calculations leading to the solution of the spatial influence 
incompressible flow for the F5 wing studied in section F of chapter m. The form of the spatial influence 
flow equation solved is introduced as Equation (3.23) and is reproduced below. The solution was obtained 
for nine equally spaced chord locations. The software utilized to perform the analytical calculations is called 
MAPLE (Ref [16]). 


1 Q( 0 ) Q(x) 1 

Q = - -^[QJ0)ln2x.QJx^)ln2(x^-x)] 

411 X X -X 47t 


471*'0 471 Jjt 

(3.23) 
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>Q:=71.2901*xi''10-352.0008*xi'"9+743.1697*xi^8-875.6007*xi''7+630.8471*xi'^6-2 
86.6580*xi''5+81.7019*xiM-14.1307*xl''3+1.4169*xi''2-0.0348*xi; 

Q ;= 71.2901 - 352.0008 + 743.1697 - 875.6007 4’ 

+ 630.8471 - 286.6580 + 81.7019 i,* - 14.1307 

+ 1.4169 ^^-.0348^ 

> clQ1:=diff(Q,xi); 

dQl := 712.9010 - 3168.0072+ 5945.3576 

- 6129.2049 + 3785.0826 - 1433.2900 

+ 326.8076 - 42.3921 ^ + 2.8338 ^ - .0348 

> clQ2:=diff(dQ1,xi); 

dQ2 ;= 6416.1090 ^ - 25344.0576 ^ + 41617.5032 
- 36775.2294 + 18925.4130 - 5733.1600 ^ 

+ 980.4228 ^ - 84.7842 4 + 2.8338 

> dQ3:=diff(dQ2,xi); 

dQ3 ;= 51328.8720 ^ - 177408.4032 + 249705.0192 

- 183876.1470 + 75701.6520 -17199.4800 

+ 1960.8456 ^ - 84.7842 

> ans:=subs(xi=:0,dQ1); 


> ans:=subs(xi=1,dQ1); 


> ans:=subs(xi=0,dQ2); 


ans := -.0348 


ans := .0536 


ans := 2.8338 
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> ans:=subs(xi=1,dQ2); 


ans := 5.0506 

> term1 :=(1/(4*Pi))*((0.0348/x)+(0.0536)/(1-x)); 

1 1 

.0348 - + .0536- ' 

_ 1 X \ - X 

term! :=- 

4 71 


> term2:=:(1/(4*PI))*(2.8338*(ln(2*x))+5.0506*ln(2*(1 -x))); 

^ 1 2.8338 ln(2;c) +5.0506 ln(2-2x) 

termz :=-—- 

4 7C 


> term3:=dQ3*ln(2*(x-xi)); 

terms := (51328.8720 - 177408.4032 + 249705.0192 

- 183876.1470 + 75701.6520 - 17199.4800 

+ 1960.8456 ^-84.7842) ln(2ji:-2^) - 

> term4:=dQ3*ln(2*(xi-x)); ~~ ' 

term4 ;= (51328.8720 ^ - 177408.4032 'C + 249705.0192 

- 183876.1470 + 75701.6520 f - 17199.4800 

+ 1960.8456^- 84.7842)ln(2 %-2x) 

>term5:=(1/(4*Pi))*lnt{term3,xi=0.0..x); 


termS := 


f Cx 

= ^ ( 

4 ^ 

*/ _ 


51328.8720 - 177408.4032 


+ 249705.0192 - 183876.1470 + 75701.6520 

- 17199.4800 ^ + 1960.8456 % - 84.7842) ln(2 jc - 2 



> term5:=(1/(4*Pi))*int(term3,xl=0.0..x); 


1 

terms := - (84.78420000 x + 83970.10713 x^ - 1470.634200x^ 
4 

+ 10510.79333 x^ - 39427.94375 x'* 

- 36775.22940 ln(2. x) x^ - 101962.8828 x* 

+ 65713.52078 x’ - 25344.05760 ln(2. x) x' 

- 84.78420000 ln(2. x) x + 41617.50320 ln(2. x) xS 

- 17438.06768 x* + 980.4228000 ln(2. x) x^ 

- 5733.160000 ln(2. x) x^ + 6416.109000 ln(2. x) x* 

+ 18925.41300 ln( 2. x)x^)/jt 


> term6:=(1/(4*Pi))*lnt(term4,xi=x..1); 

l^r. 

term6 := — 

4 


(51328.8720 - 177408.4032 


J 

V ^ 


+ 249705.0192 ^ - 183876.1470 + 75701.6520 ^ 


- 17199.4800 'C + 1960.8456 ^ - 84.7842) ln(2 ^ - 2 x) ^ 

\ 



) 

> term6:=(1/(4*Pi))*int(term4,xi=x..1); 
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1 

term6 := - (-6416.109000 ln(2. - 2. x) x* + 17438.06768 x* 

- 329.3686900 x - 41617.50320 ln(2. - 2. x) x« 

+ 3594.581720 x^ - 980.4228000 ln( 2. - 2. x) x^ 

+ 62238.41350 x^ - 115054.2845 x5 + 124098.8859 x*^ 

- 19868.30010x3- 18925.41300ln(2.-2.x)x^ 

+ 36775.22940 ln(2. - 2. x) x^ 

+ 5733.160000 ln(2. - 2. x) x^ - 72129.62978 x’ 

+ 84.78420000 ln(2. - 2. x) x + 11.63426214 
+ 2.216800000 ln(2. - 2. x) + 25344.05760 ln(2. - 2. x) x’ 

)ln 

> uoverUinf:=tenn1-term2-term5+term6; 


uoverUinf:= 


1 1 

.0348 — -I- .0536- 

^_X _ 1 -X 

4 K 


1 2.8338 1ii(2x)-h 5.0506 ln(2-2x) 1 
4 K 4 

84.78420000 X-I- 83970.10713 x^ - 1470.634200 x^ 
+ 10510.79333 x^ - 39427.94375 x* 

- 36775.22940 ln(2. x) x^ - 101962.8828 

+ 65713.52078 x’ - 25344.05760 ln(2. x) x^ 

- 84.78420000 ln(2. x) x + 41617.50320 ln(2. x) x« 
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- 17438.06768 ;c» + 980.4228000 ln( 2 . x) 

- 5733.160000 ln( 2 . x ) + 54 j ^ j 99000 ln( 2. jc) jc« 

+ 18925.41300 ln(2. ;c);c'*)/7C + ^( 

-6416.109000 ln( 2 . - 2. x) jc* + 17438.06768 x* 

- 329.3686900 X- 41617.50320 ln(2. - 2. x) x« 

+ 3594.581720x2 - 980.4228000 ln( 2. - 2. x)x2 

+ 62238.41350 x^ - 115054.2845 x^ + 124098.8859 x‘ 

- 19868.30010x2- 18925.41300 ln(2. - 2.x)x^ 

+ 36775.22940 ln( 2 . - 2 . x) x^ 

-I- 5733.160000 ln(2. - 2. x) x^ - 72129.62978 x’ 

-I- 84.78420000 ln( 2. - 2. x) x - 1 -11.63426214 

-I- 2.216800000 ln(2. - 2. x) + 25344.05760 ln(2. - 2. x) 




> xvec:=array(1.. 10); 
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> for i from 1 to 10 do uovrUinf[i]:s::evalf(subs(x=xvec[i],uoverUinf)) od; 

uovrUinf^ ^ ^ := -.0573921826 
uovrUinf^^^ := .0346121370 
uovrUinf^^^ := .0017849588 
uovrUinf := -.00548285722 

L ^ J 

uovrUinf^^^ :=-.01293665332 
uovrUinf^^^ := -.08484371084 
uovrUinf := -.1297571415 

[7 ] 

uovrUinf, 1189657140 

[ 8 ] 

uovrUinf^^^ := -.0926006053 

Error, division by zero 

> 
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APPENDIX B. TRIANGULAR WING SPATIAL INFLUENCE FLOW 
CALCULATION 

This appendix shows all the intermediate calculations leading to the solution of the spatial influence 
incompressible flow for the triangular wing studied in section G of chapter ID. Two forms of the spatial 
influence flow equations were solved and were introduced as Equation (3.23) and Equation (3.26). They 
are respectively reproduced below. Equation (3.23) is the general analytical expression, and Equation (3.26) 
is a specific Equation derived from Equation (3.23) for the triangular wing. The solutions were obtained for 
nine equally spaced chord locations. The expression for the wing cross-sectional area was provided in Ref. 
[15].The software utilized to perform the analytical calculations is called MAPLE (Ref [16]). 


1 0(0) Q(x) 1 

- ^[QJ0)ln2x.QJx^)ln2(x^-x)] 

X X ^X ^TC 

— QJA)ln2(x-l)d% . QJ,l)ln2(l-x)dl 

471 •'0 471 •'jc 

(3.23) 

^ —(l-2x) - —(l-3x)ln4x(l-x) (3.26) 

6o K 1-X 7t It 
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> Q:=4*delta*sigma*xi'^2*(1-xl); 


> dQ1:=diff(Q,xi); 


> dQ2:=diff(dQ1,xi); 


> dQ3:=diff(dQ2,xi); 


g:=4 5a4^(l -4) 


dQl :=86o|(l -^)-48aV 
dQ2 ;=86o(l-4)-16 6cy^ 


> ans:ssubs(xb1,dQ2); 


dQ3 := -24 6 o 


ans := -16 6 a 


> termi :=(1/(4*Pi))*((-4*delta*slgma)/(1-x)); 


terml := - 


k( 1 -x) 


> term2:=(1/(4*Pi))*(8*delta*sigma*(ln(2*x))-16*delta*slgma*ln(2*(1 -x))); 

1 85oln(2jc)-166aln(2-2x) 

terml :=- 

4 K 


> term3:=dQ3*ln(2*(x-xi)); 

term3 := -24 6aln(2x-2^) 


> term4:=dQ3*ln(2*(xi-x)); 


term4 := -24 5 a ln( 2 ^ - 2 jc) 


> term5:=(1/(4*Pi))*mt(term3,xi=:0.0..x); 


term5 := 


1 -24 5 a ln(2x) X + 24 5 ax 


> term6:=(1/(4*Pi))*int(term4,xl=x..1); 
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term6 := “ (-24 6 a ln(2 - 2 jc) + 24 6 a + 24 6 a ln(2 - 2 x) jc 
-245ax)/7C 

> sigma:=1/6: 

> uoverliinf:=:term1 -term2-ternn5+term6; 

-5ln(2x)--61n(2-2x) 

1 6 13 3 

uoverUinf:= - ---^- 

671(1-x)4 n 

1 -4 6 ln( 2 x) X + 4 6 X 
4 % 

1 -4 6ln(2-2x) + 45 + 46ln(2-2x)x-45x 
4 n 
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> xvec:=array{1 .. 10); 


xvec := array( 1 .. 10, [ ]) 
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> for i from 1 to 10 do uvel[i]:=evalf(subs(x=:xvec[i],uoverUinf)) od; 

MveZ ;=.2715820284 5 

wel :=.1436123850 5 

\ 

uvel := .05338583318 6 
Mve/,„ ;=-.02562370518 5 

[4] 

Mve/ :=-.1061032954 5 

L J 

-.1997561750 5 

[ 6 ] 

«veZ := -.3245121960 5 
“veZ,,, := -.5225377152 5 

Lo J 

:= -.9694453447 5 

I" J 


Error, division by zero 


>uoverUlnf:=-(1/(1-x)/Pi)+6*(1-2*x)/Pi-2*(1-3*x)*ln(4*x)*(1-x)/Pi; 

uoverUinf := 

1 ^\-2x ^ (1 - 3 x) ln(4;c) (1 - jc) 

--+ 6 -- 2 - 

(1 - x) 71 n n 


> xvec:=array(1.. 10); 

xvec := arrayC 1 .. 10, [ ]) 

> for I to 10 do xvec[i]:=l/10 od; 


xvec^ , •= — 

[1] 10 


[ 2 ] 


:■ 

[3] 


xvec^ 


xvec^ 


xvec... := 


xvec._. := 
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4 

5 


xvec^^^: 


10 


XV^C^IO] := 1 

> for i from 1 to 10 do uvel[l]:=evalf(subs(x=xvec[i],uoverUinf)) od; 

uvel,-.= 1.541706944 
[ 1 ] 

uvel := .7934866633 

I ^ i 

uvel := .3010904523 
uvel := -.1126389811 

[ 4 ] 

uvel ;= -.5263019721 
uvel^^ := -.9993975524 

[o] 

uvel := -1.608669574 
uvel^^ := -2.530129275 

L o J 

-4.572356786 

Error, division by zero 


> 
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APPENDIX C. TRIANGULAR WING CROSS-SECTIONAL FLOW 
CALCULATION 

This appendix shows all the intermediate calculations leading to the solution of the cross-sectional 
incompressible flow for the triangular wing studied in section G of chapter IQ. The form of the spatial 
influence flow equation solved is introduced as Equation (3.17) and is reproduced below. The solution was 
obtained for nine equally spaced chord locations and the process repeated for two spanwise locations (near 
the root and near the tip). For conciseness, the intermediate spanwise locations although calculated in our 
study are not presented because of their repetitive character. The software utilized to pe'rform the analytical 
calculations is called MAPLE (Ref. [16]). 


o6 2it j(x) 


P (x^) orv'k I . P (x,0-P (x^) - 
° ln| y \ - dC 

2lt 2 It Jo ^-2 

l+'T- 

s(x) 


(3.17) 
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> Q:=4*delta*sigma*x'^2*(1-x); 


Q := 4 5 o (1 - x) 


> Qbar:=Q/(sigma*delta); 


> dQbar1:=diff(Qbar,x); 


Qbar :=4x^ (1 - x) 


dQbarl :=Sx(l -x) -4jd 

> dQbar2:=:diff(dQbar1,x); 

dQbar2 := 8 - 24 x 


> dPbar2:=-8*zbar; 


> dPzetabar2:=-8*zetabar; 


dPbarl := -8 zbar 


dPzetabarl := -8 zetabar 


> termi :=(dQbar2/(2*Pi))*(ln((1/6))+ln(x)+ln(sqrt(1 -(zbar/x)'^2))); 


terml := - 
2 


(8-24x) -ln(6) + ln(x) + In 


- zbai^ 


> term2:=(dPbar2/(2*Pi))*ln(abs((1 -(zbar/x))/(1 +(zbar/x)))); 

(I zbar 


zbar In 


zbar 


terml := -4 


>term3:=(dPzetabar2-dPbar2)/(zetabar-(zbar)); 

-8 zetabar + 8 zbar 


termS := 


zetabar - zbar 


188 



>term4:s:(1/Pi)*lnt(term3,zeta=:0.0..(x))," 


-8 zetahar + 8 zhar 
zetabar - zhar 


term4 := 


> term4:=(1/Pi)*int(term3,zeta=0.0..(x)); 


term4 := -8 — 

n 


> sigma:=1/6: 


> phi:=term1-term2-term4; 


(8 - 24x) -ln(6) +ln(jc) + ln 


zbaiP- 


zbar 


zbar In 


zbar 
1 +- 

-^ + 8 - 

n n 
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> # CALCULATIONS FOR X=1/10 

> phianaly1:=subs(x=1/10,phi); 


phianalyl := — 


zbar In 


ill 

-ln(6) + ln — +ln(/T 
14 10 


-100 zhaP ') 


1-10 zhar\ 

1-1-10 zbar I 4 1 

-+- 

71 5 71 


> zvec:=array(1 .. 10); 


zvec := aiTay( 1 .. 10, [ ]) 

> for i to 10 do zvec[i]:=6*i/600 od; 


zvec ^ ^ := - 

[ 1 ] 100 




ZVeC|3j :: 


100 




zveCj5] := 
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zv^Cr,, := — 

t6] 50 


zvec 


100 


zvec 


[ 8 ]- 25 


•= - 

t9] 100 


^^^[ 10 ]* 10 


> # the division by 6 is to get the form of the answers which has no sigma 
t the denominator 

> ans:=subs(zbar=1/100,phianaly1); 


(1/6) a 


_ -ln(6)-ln(10) + ln -J99 /lOO 

14 100 ^ 


ans := 


4 1 

+- 

5 K 

> final :=:evalf(ans)/6; 


final := -.5669238108 


> ans:=subs(zbar=:1/50,phianaly1): 

> final:=evalf(ans)/6; 
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final := -.5705043033 


> ans:=:subs(zbar=:3/100,phianaly1): 

> final :=evalf(ans)/6; 

final ;= -.5766971103 

> 

> ans:=:subs(zbar=1/25,phianaly1): 

> finai:=:evalf(ans)/6; 

final :=-.5S5S9322Sl 

> 

> ans:=subs(zbar=1/20,phianaly1); 



5 K 5 - % 5 % 

> final:=:evalf(ans)/6; 

final := -.5987749621 

> 

> ans:=:subs(zbars3/50,phianaly1); 

ans := 



71 25 71 5 7C 


5 

> final:=evalf(ans)/6; 
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final :=-.6165491328 


> 

> ans:=:subs(zbars:7/100,phianaly 1); 


14 

arts :=- 

5 

4 1 


-ln(6) - In(lO) + In 


M r- 
— /51 jlOO 
100 "' 


J 


71 


5 % 


1 



In 



3 

A 


10 



17 


V 

10 

J 


n 


> final:=evaIf(ans)/6; 

final := -.6415289300 

> 

> ans:=subs(zbar=2/25,phianaly1); 



> final:=;evalf(ans)/6; 


:=-.6789331503 


> 

> ans:=subs(zbar=9/100,phianaly 1); 
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4 1 

+- 

5 71 

> final:=evalf(ans)/6; 

final := -.7453324925 


> 

> ans:=subs(zbar=1/10,phianaly1); 




4 1 

+- 

5 K 


> final:=evalf(ans)/6; 

Error, (in In) singularity encountered 
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> # CALCULATIONS FOR X=1 

> phianaly10:=subs(x=1,phi); 


phianalylO := -8 


•ln(6) + ln( 1) + ln(/l - zhar^ 


zhar In 


1 -zhar ^ 


1 + zhar 


> for I to 10 do zvec[i]:=6*ri 0/600 od: 

> 

> ans:=subs(zbar=1/10,phianaly10); 


+ 8 - 
71 


-ln(6) + In 


ans := -8 


100 


99/100 


> final:=evalf(ans)/6; 


final := 1.182733902 


>ans:=:subs(zbar=1/5,phianaly10); 


ans := -8 


-ln(6) + ln — . 24. 25 


5 7C 


> final:sevalf(ans)/6; 


final 1.176313741 
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> 

> ans:=subs(zbar=3/10,phianaly10); 



> final:=:evalf(ans)/6; 

final := 1.165463606 




> final:=evalf(ans)/6; 

final := 1.149937580 

> 

> ans:=subs(zbar=1/2,phianaly10); 



K 



> final:=evalf(ans)/6; 
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final := 1.129341166 


> ans:=subs(zbar=3/5,phianaly10); 


ans := -8 


-ln(6) + ln -J16 J25 


5 

In — 


5 % 


> final :=evaif(ans)/6; 


final := 1.103056103 


> ans:=subs(zbar=7/10,phianaly10); 


-ln(6) + ln 


ans := -8 


100 


51 /lOO 


> final:=evalf(ans)/6; 


final := 1.070082027 


> ans:=subs(zbars4/5,phianaly10); 


+ 8 - 
n 


ans := -8 


f 1 ^ 

-ln(6) + ln 

_ J 


5 7C 


1 

+ 8 - 
n 
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> finai:=evalf(ans)/6; 


final ■- 1.028648218 


> 

> ans:=subs(zbar=;9/10,phianaly10); 



> final:=:evalf(ans)/6; 


final := .9749312040 


> 

> ans:=:subs(zbars:1,phianaly10); 


^ 4n(6) + ln(0) ^ 

ans := -8-+ 4 

K 


In 


r 

0| 



2 

J 



K K 


> final:=evalf(ans)/6; 
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> h:=2*delta*(x-abs(zysigma)*(1 -x); 
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APPENDIX D. STREAMWISE VELOCITY DISTRIBUTION TEST 
ON THE RECTANGULAR WING 


This appendix is related to the discussion in section G-2 of chapter m about the velocity 
perturbation distribution in the streamwise direction on the triangular wing. It was said for this case that 
PMARC produces erroneous results. This appendix shows that for the rectangular wing the streamwise 
velocity distribution seems to present a proper behavior. The appendix is broken down as follows. 

a. Application of the slender body theory perturbation expression to corroborate the results 
subsequently obtained with PMARC, 

b. Plot showing the streamwise velocity as a function of chordwise location for various 
spanwise locations, in the case of a rectangular wing of aspect ratio 20 with a biconvex cross-section which 
thickness is tapered from root to tip. 

c. Plot showing the streamwise velocity as a function of chordwise location for various 
spanwise locations, in the case of a rectangular wing of aspect ratio 7 with a biconvex cross-section which 
thickness is tapered from root to tip.. 

d. Plot showing the streamwise velocity as a function of chordwise location for various 
spanwise locations, in the case of a rectangular wing of aspect ratio 4.5 with a biconvex cross-section which 
thickness is tapered from root to tip.. 

e. Plot showing the streamwise velocity as a function of chordwise location for various 
spanwise locations, in the case of a rectangular wing of aspect ratio 4.5 with a NACA 0012 cross-section 
which thickness is tapered from root to tip.. 

As can be observed, the PMARC results for the biconvex cross-section cases show similar behavior 
than the one predicted by the slender body theory. Furthermore, all the distributions shown are legitimate 
as the velocities (or perturbations) become smaller toward the tip, as expected. 
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>u:=((4*tau/Pi)*(1-(0.5-x)*ln(abs((1-xyx)))*Vinf)+Vinf; 

f 

T 1 - (.5 -x) In- Vinf 

u:=A —— -+ Vinf 

n 


> xvec:=:array(1..10); 

xvec := aiTay( 1 .. 10, [ ]) 

> for i to 10 do xvec[i]:=i/10 od; 


xvec^ — 
t*] 10 


xvec := 
[ 2 ] 


:= 


xvec ^^^:: 


xvec ^^^:: 


xvec^^^ 


xvec^^^:-. 
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xvec,„, := — 

[ 9 ] 10 


xv<.Cj,o,:=l 


> tau:=0.098: 

>Vinf:=11: 

> for i from 1 to 10 do uvel[i]:=evalf(subs(x=xvec[i],u)) od; 

«vgZ := 11.16623003 

uvel := 11.80172380 

uvel ;= 12.13996012 

L J 

MveZ , := 12.31690003 

[ 4 ] 

uvel := 12.37255223 

L J 

uvel, ;= 12.31690003 

LoJ 

uvel := 12.13996012 

ttveZ,;= 11.80172380 

LoJ 

Mve/ := 11.16623003 

Error, (in evalf/ln) singularity encountered 
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Slender Body Theory - Rectangular Wing 
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Vx vs X/C at Various Percent Half-Span Locations - Biconvex Set. Thickness Taper Rect. Wing, AR=20 
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Vx vs X/C 




















APPENDIX E. SPANWISE CROSS-SECTIONAL AREA CODE 


This appendix present the MATLAB computer code created to transform the chordwise section 
definition of the F5 wing into 41 spanwise cross-sectional areas. As can be seen, the program's input 
requirements are the x, y, and z coordinates of the cross-sections. The y coordinates can be taken from the 
PMARC output file, which by virtue of its construction delineate some of the geometiy parameters. Also 
required in the input are the sweep angle and span. Section B-6-a of chapter IV details how the code 
operates. 

This code could be used on other wings than the F5 providing the input data mentioned above is 
updated as well as some other parameters scattered within the code. It could become advantageous to bring 
all variables that are now scattered in the code at the beginning of the file so it is a lot easier to adapt to other 
wings. Section B-6-b of chapter IV explains why the program is presently geared more towards the F5 
wing. 
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Printed by carrier from viper 
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Printed by carrier from viper 
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Printed by carrier from viper 




231 































232 








































APPENDIX F. F5 WING SPANWISE CROSS-SECTIONAL AREAS 


This appendix provides the results of the computer code described in section A-6 of chapter IV. It 
shows the area, radius and diameter of each of the 41 spanwise cross-sections defined along the chord of the 
F5 wing. Also, a plot is included that shows the distribution of the spanwise cross-sectional areas along the 
chord. 
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Columns 22 through 28 





































APPENDIX G. TRIANGULAR WING SPANWISE CROSS-SECTIONAL 
AREAS 


This appendix provides the radii of 46 spanwise cross-sectional areas with the x coordinate of their 
centroid along the chord of the triangular wing. This information is provided for the two cases of maximum 
thickness ratio treated: 0.098, and 0.196. A plot showing the distribution of the spanwise cross-sectional 
areas along the chord is also included for each case. 
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Spanwise Cross-sectional Area 


Triangular Wing AR=!2/3, delta=0.196, Q vs x 
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APPENDIX H. AUTOMATED EXTRACTION OF CL AND CM 


This appendix presents two short programs written in the UNIX AWK language used to extract CL 
and CM at each time-step of the PMARC output file 'DATA6'. These programs introduced in section B-4 
of chapter V are specific to the case of the F5 wing treated in this study. The location of CL and CM in the 
output file can vary depending on the configuration modelled. 
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CM EXTRACTION 


#!/bin/sh 

# AWK Script to extract the total cm value from 

# the PMARC output DATA6 file. 

# Finds "WIND AXES" and skips 32 lines then picks the fourth field ($4 

# which is the total cm value 

# 

echo 

echo Reads in DATA6 file for a simple wing with symmetry 

echo and its wingtips for a total of four patches 

echo 

echo cm total data is put in file \’'outcm\* 
awk ' 

/WIND AXES/ { 

NL = 32 

for{i=0;i<NL;++i) 

getline # Gets 32 line after WIND AXES term 

print $4 # Gets fourth field (cm value) 

} 

' DATA6 > outcm 


CL EXTRACTION 

***★★*•★★★*★★■* ^ 

#i/bin/sh 

# AWK Script to extract the total CL value from 

# the PMARC output DATA6 file. 

# Finds "WIND AXES" and skips 32 lines then picks first field ($1) 

# which is the total CL value 

# 

echo 

echo Reads in DATA6 file for a simple wing with symmetry 

echo and its wingtips for a total of four patches 

echo 

echo CL total data is put in file \"outCL\* 
awk ' 

/WIND AXES/ { 

NL = 32 

for(i=0;i<NL;++i) 

getline # Gets 32 line after WIND AXES term 

print $1 # Gets first field (CL value) 

} 

' DATA6 > outCL 
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